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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 


The  purpose  of  the  work  reported  here  was  to  investigate  nonsteady-state  behavior 
in  the  early  stages  of  one-dimensional  shock  propagation.  This  behavior  has  been 
termed  stress  relaxation.  The  work  consisted  of  an  experimental  phase  and  a 
theoretical-computational  phase.  In  the  experimental  pha^e,  stress-time  profiles 
resulting  from  impact  by  a  gas-gun  projectile  were  observed  in  hard  and  soft  2024 
aluminum  and  in  polyethylene  at  various  distances  from  the  ^lane  of  impact. 

Quartz  gages  were  used  for  the  measurements  in  aluminum,  and  manganin  wire  gages 
were  used  in  polyethylene.  Peak  stress  amplitudes  were  about  14  kbar  in  aluminum; 
there  was  a  pronounced  decay  in  amplitude  of  elastic  precursoi  with  propagation 
distance.  The  peak  stress  in  polyethylene  was  varied  from  3  to  11  kbar;  non- 
steady-state  effects  were  also  observed,  although  no  elastic  wave  was  evident. 

In  the  theoretical  and  computational  phase,  several  mathematical  models  of  stress 
relaxation  were  developed.  Computations  based  on  these  models  have  been  performed 
by  »~he  method  of  characteristics  and  by  the  artificial  viscosity  method  of  von 
Neumann  and  Richtmyer.  It  has  been  found  that  the  two-parameter  model  of  stress 
rel *  .ation  predicts  the  experimental  results  more  closely  than  does  the  one- 
parameter  model.  During  this  investigation,  it  was  found  that  no  artificial 
viscosity  is  needed  when  stress-relaxing  models  are  used.  Some  preliminary 
results  on  the  inclusion  of  thermodynamic  behavior  have  been  obtained,  and  it 
is  shown  that  a  consistent  model  can  be  formulated  on  the  basis  of  equilibrium 
thermodynamics . 
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SECTION  I 


INTRODUCTION 


:sf 

One  of  t!ie  current  tasks  of  the  Air  Force  Weapons  Laboratory  'IaFWL) 
is  to  achieve  a  physical  understanding  and  mathematical  description  of 
shock  wave  propagation  in  solids.  This  is  a  major  undertaking  because  of 
the  breadth  of  interest  in  materials,  pressure  ranges,  and  time  scales. 
Stanford  Research  Institute  has  attempted  for  several  years  to  aid  the  Weap¬ 
ons  Laboratory  in  this  task  by  carrying  out  research  on  shock  propagation, 


equations  of  state ,  shock  measurement  techniques,  and  structural  response.  The 
work  reported  here  is  a  continuation  of  our  eifort  to  be  responsive  to  the 
needs  of AFWL.  We  are  concerned  here  with  the  early  stages  of  shock  formation 
and  in  particular  with  the  phenomenon  of  “stress  relaxation.”  This  term 
refers,  in  the  present  context,  to  a  change  with  time  of  the  shear  yield 
strength  of  a  shock- compressed  material.  The  time  dependence  of  the  shear 
strength  arises  from  the  combination  of  the  very  rapid  loading  (strain) 
rate  and  the  finite  time  required  for  the  material  to  flow  under  the  ap¬ 
plied  shear  stress. 


To  describe  shock  wave  propagation,  the  equations  of  motion  must  be 
supplemented  by  constitutive  relations  describing  the  physical  state  of 
the  material,  e-g.,  the  s t r ess- s t r ai n  behavior.  The  constitutive  relations 
in  general  use  at  present  are  considered  to  be  independent  of  the  time 
scale  of  the  problem.  This  treatment  can  be  only  an  approximation,  however, 
because  real  solids  require  a  finite  time  to  respond  to  an/  stimulus  and 
hence  will  show  time-dependent  properties  on  some  time  scale. 

In  tnis  project  all  the  experiments  were  performed  using  plane  geom¬ 
etry,  and  the  computational  models  were  set  up  to  describe  plane  wave 
propagation.  Let  us,  therefore,  consider  the  plane  case  in  some  detail. 

The  initial  shock  front  is  produced  by  a  very  rapid  acceleration  of  a 
.surface  of  the  specimen  (in  the  x  direction,  say),  as  by  impact  of  a  pro¬ 
jectile.  If  the  initial  acceleration  is  infinite  then  one  may  consider 
the  initial  compression  (along  ABC  in  Figure  1,  for  example)  of  the  material 
to  he  elastic..  Because  the  strain  must  be  uniaxial,  i.e.,  the  macroscopic 
strain  is  zero  parallel  to  the  shock  front,  and  because  the  material  can 


1 


jpport  a  shear  stress,  the  lateral 

stress  cj  (=cr  )  will  be  less  then  the 
y  ‘ 

normal  stress  cr  .  However,  the  in- 
siantemeous  shear  stress  produced 
is  greater  than  the  elastic  forces 
can  support  and  shear  yielding  begins 
on  a  microscopic  scale — even  though 
macroscopic  lateral  strain  remains 
zero.  This  yielding,  requiring  a 
finite  time,  allows  a  to  relax  from 
its  instantaneous  value  at  C  to  an 
equilibrium  value  at  D  as  the  shear 
stress  relaxes  to  its  equilibrium 
value.  The  state  D  is  now  reached 
a  ross  two  waves,  the  first  being 
an  elastic  precursor  wave  of  final 
amplitude  AB  (B  is  called  the 
Hugoniot  elastic  limit),  and  the 
second  having  amplitude  BD.  During 
the  relaxation  process,  the  precursor  amplitude  will  decrease  with  propaga¬ 
tion  distance,  asymptotically  approaching  AB.  The  proximity  of  point  D  to 
the  curve  of  hydrostatic  compression  (isothermal)  depends  on  the  magnitude 
of  the  equilibrium  shear  strength.  Tf  the  material  can  support  essent’ally 
zero  shear  stress  after  reaching  equilibrium,  point  D  lies  on  the  hydrostat. 
The  larger  the  equilibrium  shear  strength,  the  higher  point  D  lies  above 
the  hydrostat.  Equilibrium  stresses  greater  than  those  predicted  from 
quasi-static  tensile  tests  are  frequently  observed.  This  phenomenon  may 
occur  because  the  eq-'i 1 i hri urn  yield  strength  depends  upon  the  hydrostatic 
stress.  On  the  other  hand  it  may  be  only  an  apparent  equilibrium  value 
because  the  relaxation  time  is  long  (many  microseconds)  compared  with  the 
duration  of  a  shock  experiment. 

Section  II  of  this  report  describes  exploratory  stress-relaxation 
experiments  in  aluminum  and  polyethylene.  This  work  was  done  in  the 
Poulter  laboratories  of  Stanford  Research  Institute.  Section  III  is  con¬ 
cerned  with  the  development  of  mathematical  models  and  techniques  to 
describe  stress  relaxation.  This  work  was  performed  in  the  Institute's 
Mathematics  Department.  Although  the  two  phases  were  carried  out  rela¬ 
tively  independently,  the  experimental  phase  defined  the  phenomena  to  be 


FIG.  1  ILLUSTRATION  OF  RELAXATION 
FROM  INSTANTANEOUS  STRESS 
(Point  C)  TO  EQUILIBRIUM  STRESS 
(Pcint  D) 


described  mathematically  and  also  provided  s  Tie  data  to  be  used  in  ap¬ 
plying  the  theory  that  was  developed. 

This  report  builds  on  previous  (nonsystematic )  studies  in  these  and 
other  laboratories;  however,  the  problems  investigated  are  in  no  sense 
completely  solved.  Much  more  work  will  be  required  before  time-dependent 
effects  in  shock  propagation  are  fully  understood. 


SECTION  II 

EXPERIMENTAL  WORK 


by 

G.  D.  Anderson  and  W.  J.  Murri 


1.  Background 

The  purpose  of  the  experimental  phase  of  this  program  was  to  look 
for  nonsteady  state  effects  in  the  early  stages  of  shock  propagation  and 
to  provide  data  as  a  guide  in  the  choice  of  models  to  be  used  in  the 
computational  phase.  Two  materials,  one  metal  and  one  plastic,  which 
would  be  expected  to  exhibit  different  types  of  nonsteady  behavior  were 
chosen  for  study.  The  metal  was  2024  aluminum.  Since  the  Hugoniot 
elastic  limit  of  2024  aluminum  is  known  to  depend  upon  the  hardness, 
both  very  soft  and  very  hard  samples  were  studied.  The  plastic  chosen 
was  a  highly  crystalline  polyethylene. 

For  both  materials  the  experiments  were  designed  to  observe  the 
shape  of  an  effectively  infinite  duration  plane  pressure  pulse  as  a 
function  of  distance  from  the  impact  face.  The  pulses  were  generated  by 
impact  from  a  thick,  flat  projectile  accelerated  by  a  gas  gun.  In  most 
cases  the  projectile  head  was  of  the  same  material  as  the  target  sample; 
a  constant  particle  velocity  boundary  condition  could  thus  be  maintained 
at  the  impact  surface  plane.  Stress-time  profiles  were  recorded  by 
quartz  gages  for  the  aluminum  and  by  mangenin  wire  gages  for  the 
polyethylene. 

2.  Experimental  Technique 

a .  Aluminum 

Previous  work  on  2024  aluminum  (Ref.  1)  has  demonstrated  an  elastic 
precursor  amplitude  varying  from  about  1  to  6  kbar  depending  upon  the 
initial  hardness  of  the  material.  In  this  project  we  investigated 
uonsteady-state  elastic-plastic  behavior  in  hard  and  soft  aluminum.  The 
2024  alloy  was  chosen  instead  of  pure  aluminum  since  a  much  greater 
difference  in  hardness  could  be  obtained  by  heat  treatment  of  the  alloy. 
Sheets  of  stock  2024-T351  aluminum  ranging  in  thickness  from  0.0195  to 


0.627  inch  were  commercially  heat  treated.*  The  greateat  hardness  was 
attained  by  artificial  aging  for  nine  hours  at  a  temperature  of  375°F. 

The  softest  aluminum  was  attained  by  annealing.  This  process  consisted 
of  heating  the  material  to  750°F  at  a  rate  of  50°F/hour,  holding  for  two 
hours  at  750°F,  cooling  at  50°F/hour  to  450°F,  and  then  air  cooling.  The 
hardness  of  the  aged  aluminum  was  Rockwell  B85  and  that  of  the  annealed 
aluminum  was  Rockwell  BO. 

After  heat  tieatment  the  specimens  were  cut  into  circular  plates  for 
mounting  on  target  assemblies  and  gas  gun  projectile  heads  Target  and 
projectile  flatness  are  essential  for  achieving  simultaneous  impact. 
Therefore  all  projectile  heads  and  target  plates  were  lapped;  departure 
from  flatness  was  held  to  less  than  0.0005  inch  across  a  2%-inth  diameter. 

Quartz  gages  were  chosen  to  record  the  stress-time  profiles  in 
aluminum  because  of  the  close  impedance  match  of  aluminum  and  quartz.  The 
target  assembly  construction  showing  the  sample,  qrartz  gage,  tilt  pins, 
and  epoxy  support  is  3hown  schematically  in  Figure  2.  The  four  tilt  pins 
are  positioned  within  0.002  mm  of  the  impact  surface  and  record  the 
simultaneity  of  impact.  The  epoxy  support  is  made  of  “C-7”  epoxy  to 
which  glase  beads  nave  been  added  to  increase  its  shock  impedance  to 
match  that  of  quartz,  thereby  minimizing  the  effect  of  lateral  rare¬ 
factions. 

A  complete  description  of  the  technique  of  mounting  and  aligning 
the  target  perpendicular  to  the  axis  of  the  gun  barrel  has  been  given  by 
Linde  and  Schmidt  (Ref.  2).  Aluminum  projectile  heads  of  the  same 
material  as  the  target  were  bonded  onto  stock  aluminum  projectiles  with 
"910”  cement.  Care  was  taken  to  insure  that  the  plane  of  the  flat  face 
of  the  head  was  perpendicular  to  the  axis  of  the  projectile.  The  pro¬ 
jectile  heads  were  made  sufficiently  thick  so  that  no  reflected  wave  from 
the  back  surface  of  the  projectile  head  could  reach  the  quartz  gage 
during  the  time  the  pressure  profile  was  being  recorded. 

The  projectile  velocity  was  maintained  at  approximately  600  ft/sec, 
a  speed  which  produced  a  plastic  wave  amplitude  of  about  14  kbar.  This 
choice  was  made  so  that  the  steady-state  elastic  wave  amplitude  would  be 
approximately  half  that  of  the  plastic  wave  amplitude  in  the  hardened 
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aluminum.  In  this  way  both  waves  could  be  conveniently  observed  on  the 
same  oscilloscope  trace  and  the  linear  range  of  the  quartz  crystal  gage 
would  not  be  exceeded. 

b.  Polyethylene 

All  the  polyethylene  samples  were  cut  from  the  same  1^-inch  thick 
sheet  of  “Hi-fax." *  Test  cores  were  taken  at  several  positions  in  the 
sheet  and  examined  for  uniformity  of  density  and  crystallinity.  The 
density  as  determined  by  volume  and  mass  measurements  was  0. 938  ±  0. 002  gm/cm 
X-ray  patterns  were  taken  in  three  mutually  perpendicular  directions.  The 
patterns  indicated  a  degree  of  crystallinity  greater  than  90  percent  and 
very  small  differences  in  orientation  between  directions.  For  shot  con¬ 
struction,  polyethylene  samples  were  cut  into  circular  disks  ranging  in 
thickness  from  0.031  tc  0,75  inch. 

Manganin  wire  pressure  transducers  (Ref.  3)  were  chosen  to  monitor 
the  stress  profiles  in  the  polyethylene.  The  reasons  for  this  choice 
were  that  (1)  the  impedance  mismatch  between  quartz  and  polyethylene 
would  complicate  interpretation  of  quartz  records,  and  (2)  manganin  could 
be  imbedded  directly  in  polyethylene,  thereby  introducing  a  minimum  per¬ 
turbation  on  the  wave  profile,  Manganin  wire  gages  also  occupy  much  less 
area  on  a  specimen  for  a  given  length  of  recording  time  than  do  quartz 
gages.  Because  of  their  small  size  it  was  possible  to  use  two  wires  at 
different  depths  in  the  target  from  the  impact  plane  in  each  shot.  This 
made  it  possible  to  observe  the  same  wave  front  at  two  positions. 

A  shot  assembly  is  shown  in  Figure  3.  The  usual  target  alignment 
procedure  (see  Ref.  2)  relies  upon  the  electrical  conductivity  of  the 
target;  since  polyethylene  is  nonconductive  this  procedure  was  modified. 

An  additional  O-ring  groove  was  cut  in  the  adapter  ring  so  that  the  target 
assembly  could  be  held  by  vacuum  seal  to  the  adapter  ring.  A  lapped 
aluminum  plate  was  substituted  for  the  polyethylene  assembly  and  the 
adapter  ring  was  aligned  according  to  the  procedure  outlined  in  Ref.  2. 

The  adapter  ring  was  then  clamped  securely  into  position  to  prevent,  move¬ 
ment  during  venting  of  the  gun  barrel  to  almoapheric  pressure.  The  flat 
aluminum  plate  was  then  removed  and  replaced  with  the  polyethylene  target 
assembly.  Finally,  the  gun  barrel  was  evacuated  again. 
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Polyethylene  heads  were  mounted  on  aluminum  projectiles.  The  poly* 
ethylene  heads  were  0.75  inch  thick  to  prevent  reflected  waves  from  the 
aluminum  backing  from  overtaking  the  initial  shock  during  the  experiment. 
In  order  to  clo~*  the  tilt  pins  upon  impact  it  was  necessary  to  cover  the 
polyethylene  projectile  head  with  a  thin  layer  (0.001  inch)  of  aluminum 
foil  which  was  electrically  connected  v.o  the  projectile  body  and  hence  to 
ground.  Flat  surfaces  were  much  more  difficult  to  prepare  on  polyethylene 
projectile  heeds  than  they  were  on  aluminum.  Consequently  both  projectile 
tilt  and  nonflatness  of  the  target  contributed  to  nonsimultaneity  of 
impact.  The  contribution  to  wave  tilt  due  to  projectile  tilt  was  reduced 
significantly  by  using  12-inch  long  projectiles  weighing  approximately 
2400  gm,  whereas  3-inch  projectiles  weighing  about  260  gm  were  used  in  the 
aluminum  studies. 

The  peak  pressures  induced  into  most  of  the  polyethylene  targets  were 
in  the  2  to  3  kbar  range,  although  three  shots  were  performed  at  higher 
pressures . 


c.  Recording  Voltage  Acrosa  Manganin  Wire  at  Low  Stress  Levels 

The  voltage  change  across  the  sensitive  portion  of  the  manganin  wire 
for  shock  pressures  of  2  to  3  kbar  is  10  to  20  mv .  It  is  necessary  to 
observe  this  voltage  variation  on  top  of  a  2v  turn-on  step  in  rhe  sensing 
wire.  As  discussed  in  Ref.  3,  the  ratio  of  the  change  in  resistance  of 
the  wire  to  the  initial  resistance  of  the  wire  is  proportional  to  the 
stress. 

A  constant  current  (turned  on  a  few  si c roseconda  prior  to  shock 
arrival)  is  passed  through  the  wire  during  an  eAp^riment  and  the  percent¬ 
age  change  in  voltage  drop  across  the  sensitive  portion  of  the  wire  is 
measured.  Consequer  ly  it  is  necessary  to  measure  either  the  vslue  of 
the  current  or  the  initial  voltsge  across  the  wire  (approximstely  2v), 
and  to  measure  the  change  in  voltsge  due  to  the  passage  of  the  shock  wave 
(approximately  10  to  20  mv).  The  pressure  as  a  function  of  time  can  be 
found  from  the  resulting  oscilloscope  records. 

In  practice  certain  oscilloscopes  have  their  sensitivity  adjusted 
vo  observe  the  entire  voltage  change  across  the  wire,  the  2v  turn-on 

plus  the  10-20  mv  pressure  signal  on  top  of  it;  other  oaci  lloscopes 
operating  at  greater  aenaitivity  observe  only  the  10-20  mv  preaaure 
voltage  in  order  to  resolve  the  detail  of  the  wave  sh/ipe.  To  display  the 
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signal  due  to  the  small  pressure  pulse,  a  negative  2v  must  be  fed  into 
the  oscilloscope  to  balance  the  positive  2v  from  the  power  supply.  The 
power  supply  voltage,  however,  is  not  always  exactly  2v  (or  whatever  the 
predetermined  value  is),  but  may  vary  by  a  few  percent  depending  upon 
the  length  of  time  the  power  supply  capacitors  are  allowed  to  charge. 

This  small  variation  between  the  power  supply  turn-on  voltage  and  the 
preset  compensating  voltage  may  be  sufficient  to  cause  the  sweep  of  the 
high-sensitivity  oscilloscopes  to  be  off-scale.  Therefore  a  special  cir¬ 
cuit  was  designed  to  eliminate  the  need  for  presetting  a  compensating 
voltage  to  offset  the  power  supply  voltage. 

The  circuit  is  shown  schematically  in  Figure  4.  Any  signal  from  the 
gage,  such  as  the  turn-on  signal  from  the  power  supply,  is  sent  along  two 


FIG.  4  TIME-DELAY  REFLECTOMETER  CIRCUIT 


parallel  paths.  On  one  path  it  goes  immediately  to  the  A  input  of  the 
amplifier  and  on  to  the  oscilloscope.  Along  the  other  path  it  is  delayed 
for  2  usee  and  is  then  fed  into  the  negative  B  input  of  the  amplifier, 

At  2  c  after  turn  on,  the  circuit  returns  the  signal  reaching  the  os¬ 
cilloscope  to  sero.  Thus  if  the  power  supply  is  turned  on  at  least  2  u-**c 
prior  to  the  pressure  signal,  the  turn-on  voltage  will  be  autoaiatical 1 / 
compensated  for,  regardless  of  its  amplitude.  The  circuit  then  provides  a 
time  window  of  2  usee  duration  for  observing  tbe  pressure  pulse  signal  be¬ 
fore  it  too  is  cancelled.  The  delay  time  may  be  varied  by  varying  the 
length  of  cable  in  the  delay  line. 

3.  Experimental  Results  and  Conclusions 

Ali  shots  fired  during  the  project  are  listed  by  number  snd  objective 
in  Appendix  1.  Shots  which  were  of  an  exploratory  or  development  nature 
are  not  further  discussed  in  the  following  presentation  of  data. 


a.  A1 uminum 


The  experimental  data  for  hardened  and  annealed  aluminum  are  pre¬ 
sented  in  Table  I.  Shock  velocities  computed  from  the  time  interval 
between  the  firing  of  the  tilt  pins,  which  denotes  impact,  and  the  ap¬ 
pearance  of  a  portion  of  the  wave  at  the  quartz  gage  are  generally  lower 
than  the  accepted  steady-state  values  for  aluminum.  The  discrepancies 
appear  greater  for  the  thinner  samples.  Elastic  velocities  are  deter¬ 
mined  from  the  time  interval  between  projectile  impact  and  arrival  of  the 
first  disturbance  at  the  quartz  gage.  A  0. 0001-inch  gap  between  the 
quartz  gage  and  the  aluminum  specimen  resulting  from  a  deviation  in  spec¬ 
imen  flatness  could  produce  a  delay  in  the  response  of  the  quartz  gage 
which  would  result  ir.  the  observed  discrepancies.  The  major  purpose  of 
the  experiments  was  to  observe  wave  shape.  Had  wave  velocity  been  of 
primary  interest  different  measurement  techniques  would  have  been  used. 

The  stress  wave  amplitudes  in  the  aluminum  were  taken  to  be  those 
recorded  by  the  quartz.  No  Hugoniot  data  in  the  10  to  20  kbar  range  for 
aluminum  of  the  type  used  here  are  available  to  form  a  basis  for  estimating 
the  quart r- al umi num  mismatch.  Graphical  interpolation  of  aluminum  data 
(Ref.  1)  between  the  Hugoniot  elastic  limit  and  about  25  kbar  indicates  no 
difference  in  the  impedance  of  quartz  and  that  of  soft  aluminum.  The 
hardened  aluminum  i  slightly  higher  in  impedance;  this  could  result  in 
the  quartz  reading  as  much  as  0.5  kbar  lower  than  the  true  aluminum  pres¬ 
sure  at  the  15  kbar  level.  Since  the  effects  observed  a:e  t ime- dependent , 
the  st r ess- part ic 1 e  velocity  state  in  the  aluminum  does  not  lie  op  the 
Hugoniot.  Therefore,  since  the  impedance  mismatch  is  very  small,  an  at¬ 
tempt  to  correct  the  quartz  reading  using  an  interpolated  graphical 
Hugoniot  for  aluminum  when  the  aluminum  state  does  not  lie  on  that  Hugoniot 
would  not  significantly  enhance  the  value  of  the  data. 

Have  profiles  in  hardened  aluatinum  as  recorded  by  quartz  gages  are 
shown  in  Figure  5.  The  sample  thicknesses  varied  from  3.2  to  15,7  mm. 

In  Figure  5a.  the  record  from  a  3.2  bub  sample,  it  can  be  seen  that  a  pre¬ 
cursor  whose  amplitude  was  nearly  that  of  the  plastic  wave  was  forming. 

In  Figures  5b  and  5c  (records  from  7.89  mm  and  15-75  mm  samples),  the  pre¬ 
cursors  have  separated  more  from  the  plastic  wave  and  decreased  in  aa^ili- 
tude.  The  Hugoniot  elastic  limit  for  2024-T4  is  about  5.4  kbar.  The 
aluaiinum  studied  here  is  considerably  harder  than  T4,  and  the  quartz  record 
indicates  a  precursor  amplitude  of  7.9  kbar  at  a  depth  of  15.15  mm.  It  is 
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Table  I 

ALUMINUM  SHOT  DATA 


SHOT 

NUMBER 

PROJECTILE 

TARGET 

VELOCITY 
OF  FIRST 
MOTION 

ELASTIC  *AVEk 

,  PLASTIC  WAVE6 

Tilt* 

(^sec ) 

Velocity 

(■./“■eel 

Type  of 
A1  uninun 

Thickness 

(na) 

A 

(•a/Mec) 

A 

(m/Mcc) 

A 

(kbar) 

A 

(aa/Micc ) 

A 

(kbar) 

A 

(kbar) 

12172 

0.032 

0.171 

Soft 

3.22 

5.94 

- 

- 

- 

5.39 

- 

12.6 

12168 

0.025 

0.197 

Soft 

7.91 

5.76 

5.34 

1.6 

2.9 

5.0 

- 

15.5 

12170 

0.043 

C.  181 

Soft 

7.92 

6.00 

- 

- 

- 

5.45 

- 

13.0 

12175 

0.016 

0.205 

Soft 

7.93 

6.23 

5.90 

2.3 

3.8 

5.58 

- 

14.6 

12169 

0.060 

0.194 

Soft 

15.86 

6.00 

5.83 

1.4 

2.0 

~5.0 

- 

12.5 

12807 

0.014 

0.197 

Herd 

0.475 

- 

- 

- 

- 

- 

- 

14.0e 

12808 

0.043 

0. 197 

Hard 

0.788 

- 

- 

- 

- 

- 

- 

14.  ^ 

12173 

<0.010 

0.  184 

Hard 

3.26 

5.3 

4.27 

12.2 

12.6 

'■'■•3.5 

- 

13.1 

12171 

0.048 

0.180 

Hard 

7.89 

6.02 

5.62 

10.0 

10.7 

5.0 

- 

13.5 

12174 

0.015 

0.202 

Hard 

15.85 

6.29 

6.10 

7.1 

7.9 

5.41 

14.6 

16.1 

*  Tilt  icrosi  center  conductor  of  quarts  gage  (H  inch). 

**  See  diagram  for  definition  of  symbols. 

c  Strea*  relaxed  to  13.0  kbar  0.5  y,*ec  after  inpact  (0.31  u«ec  after  peak  »tret»s  ess  reached). 
St reaa  subsequently  increased  beyond  1  4.0  kbar. 

Stress  relaxed  to  13.2  kbar  0.t>u**c  after  inpsct  (0.35  usee  after  peak  street  ess  reached). 
Stress  subsequently  increased  beyond  ’  4. 2  kbar. 


(a)  SHOT  12,173  - 3.261  mm  THICK 


(b)  SHOT  12.i7l - 7.899mm  THICK 


i - • 
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(C)  SHOT  12.174 - 15.85mm  THIO* 


FIG.  5  QUARTZ  GAGE  RECORDS 

FOR  HARDENED  2024  ALUMINUM 


believed  that  at  this  depth  the  elastic  wave  amplitude  had  closely  ap¬ 
proached  the  Hugoniot  elastic  limit. 

Wave  profiles  in  soft  2024  aluminum  show  a  different  type  of  behavior. 
Quartz  gage  records  are  shown  in  Figure  6.  As  with  the  hard  aluminum,  the 
projectile  velocity  was  maintained  at  about  600  ft/sec,  inducing  a  plastic 
»  wave  amplitude  of  about  14  kbar  into  the  material.  The  quartz  record  in 

Figure  6a  is  the  wave  profile  emerging  from  a  3.2  mm  sample.  A  dip  in  gage 
voltage  occurred  after  the  initial  rise  to  near-peak  plastic  wave  pressure. 
This  dip  has  been  observed  on  several  shots  (see  the  elastic  wave  profile 
in  Figure  5c).  It  is  not  presently  being  interpreted  as  a  precursor  be¬ 
cause  its  duration  and  amplitude  were  always  approximately  the  same  (where- 
ever  it  occurred),  regardless  of  specimen  thickness  and  hardness. 

This  dip  at  top  of  the  quartz  stress  profile  has  also  been  seen  witn 
quartz  on  “ Plexiglas, ”  so  it  is  not  clear  at  present  whether  this  is  due  to 
the  behavior  of  quartz  under  certain  circumstances  or  whether  it  is  a 
nonsteady- state  effect  in  the  materials  under  consideration.  Excluding 
the  dip  at  r.he  top  of  the  record  in  Figure  6a,  r.o  indication  of  a  double 
wave  is  evident.  The  quartz  record  from  a  7.94  mm  thick  specimen  is  shown 
in  Figure  6b.  Here  again  the  dip  in  the  gage  record  appears  near  the  top 
of  the  wave  profile.  However,  a  precursor  is  clearly  evident  near  the 
bottom  of  the  pressure  rise.  In  Figure  6c,  the  quartz  record  shows  the 
stress  profile  emerging  from  a  15.86  mm  sample;  the  peak  pressure  is  not 
•  visible  but  the  precursor  is  clearly  well  out  in  advance  of  che  plastic  wave. 

The  wave  profiles  ir.  the  soft  aluminum  differ  from  those  in  the  hard 
in  that,  when  the  precursor  is  resolvable  from  the  plastic  wave,  it  has 
slmost  achieved  its  steady- st»te  amplitude.  The  precursor  in  the  hard 
aluminum,  however,  breaks  out  at  an  amplitude  only  slightly  less  than  that 
of  the  plastic  wave.  This  implies  that  the  relaxation  of  elastic  stress 
occui.'  at  a  more  rapid  rate  in  soft  annealed  aluminum  than  it  does  in 
hardened  aluminum.  Relaxation  of  shear  stress  is  achieved  by  plastic  flow. 

It  is  quite  reasonable  to  expect  thi s  flow  to  occur  more  rapidly  in  annealed 
aluminum,  where  the  dislocations  are  more  mobile,  than  in  hardened  aluminum 
where  the  dislocations  are  pinned.  The  decay  of  precursor  amplitude  with 
propagation  distance  ia  shown  in  Figure  7. 

In  a  symmetrical  impact  of  a  thick  projectile  into  a  thick  target  of 
s'  the  same  material,  a  relaxation  of  the  plasti*'  wave  amplitude  might  also 

)j  1  be  expected.  Consider  a  material,  such  as  aluminum,  which  behaves 
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{ a  }  SHOT  12,172 —  3.23mm  THICK 


( b )  SHOT  12,175  —  7.94mm  THICK 


( c  )  SHOT  12,169  —  1586mm  THICK 


FIG.  6  QUARTZ  GAGE  RECORDS  FOR 

FOR  ANNEALED  2024  ALUMINUM 
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FIG.  7  DECAY  OF  ELASTIC  PRECURSOR  AMPLITUDE 

WITH  PROPAGATION  DISTANCE  IN  2024  ALUMINUM 


elastically  up  to  a  certain  normal  stress  cr t  in  uniaxial  strain  and  then 
yields  by  plastic  flow.  The  stress-particle  velocity  Hugoniot  is  shown 
in  Figure  8  bv  the  curve  OBC,  and  the  stress  O’ ,  above  which  the  material 
no  longer  behaves  elastically  under  steady-state  conditions,  is  called 
the  Hugoniot  elastic  limit.  If  this  material  is  impacted  by  a  projectile 
of  the  same  material  moving  at  a  velocity  If,  the  stress-particle  velocity 
state  in  the  projectile  when  steady  state  is  reached  will  be  on  the  curve 
WB'C ■  Because  of  the  symmetry  of  the  impact  the  curves  OBC  and  WB'C  are 
mirror  images  reflected  about  particle  velocity  W/2. 

Let  us  now  assume  that  the  impact  is  perfectly  simultaneous  (no  tilt) 
and  that  the  response  of  a  mass  element  to  a  step  input  is  initially 
elastic  followed  by  a  relaxation  of  the  excess  shear  stress  above  the 
yield  stress  to  the  final  plastic  state.  Lpon  impact  the  material  at  the 
interface  is  initially  shocked  to  state  A  bicause  of  its  elastic  response. 
As  relaxation  of  the  shear  stress  occurs  the  stress  decreases  and  roaches 
C  when  relaxation  is  complete.  The  relaxation  of  the  stress  at  the  inter¬ 
face  is  communicated  to  the  wave  fro.it,  which  during  this  relaxation  time 
will  have  propagated  into  the  material.  The  " attenuation "  of  peak  initial 
stress  due  to  the  relaxation  results  in  some  interior  particle  being 
shocked  elastically  to  a  lower  amplitude  stress  state,  such  as  D.  This 
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PARTICLE  VELOCITY 

FiG.  8  PLOT  OF  STRESS  v.  PARTICLE  VELOCITY  SHOWING  A  POSSIBLE 
MECHANISM  FOR  STRESS  RELAXATION 

particle  then  relaxes  toward  C  along  some  path.  All  particles  must  ulti¬ 
mately  reach  state  C  because  of  the  symmetry  of  the  impact.  As  yet  no 
double  wave  structure  has  been  formed.  At  some  point  the  stress  at  '.he 
wave  front  will  have  been  attenuated  to  point  E,  which  is  just  the  value 
of  the  steady-state  final  stress.  The  mass  element  which  is  initially 
shocked  to  E  achieves  sts  relaxation  by  accelerating  at  constant  stress 
to  point  C.  As  a  result  of  the  stress  and  particle  velocity  gradients  in 
the  wave,  some  subsequent  particle  is  shocked  elastically  to  point  F  and 
is  then  further  compressed  plastically  to  state  C.  It  is  at  states  below 
E  that  the  double  wave  structure  appears.  The  particle  that  is  shocked 
through  states  GFC  experiences  a  double  wave  structure.  This  process  of 
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decay  of  the  elastic  precursor  continues  cntil  point  B,  the  Hugoniot 
elastic  limit,  is  reached.  Particles  at  greater  depth  experience  two 
steady-state  waves  in  which  they  are  shocked  elastically  to  B  and  then 
plastically  to  C. 

This  is  just  one  suggested  model  for  the  process  by  which  the  wave 
could  approach  steady  state.  Other  equally  plausible  models  can  be  con¬ 
sidered,  in  which  the  initial  states  lie  on  curves  other  than  linear 
extensions  of  the  elastic  portion  of  the  Hugoniot.  However,  in  all 
models  the  boundary  condition  of  constant  particle  velocity  at  the  inter¬ 
face  is  required  by  symmetry  of  the  impact.  This  means  that  the  inter¬ 
face  relaxes  from  the  initial  stress  to  point  C  along  a  path  of  constant 
particle  velocity.  Other  particles  in  the  target  will  accelerate  with 
decreasing  or  increasing  stress  toward  point  C. 

In  the  experiments  no  significant  decrease  in  the  plastic  stress 
was  observed.  The  final  plastic  stresses  recorded  were  consistent  with 
the  stress  predicted  from  existing  Hugoniot  data  and  the  measured  pro¬ 
jectile  velocities.  An  -"xtension  of  the  elastic  portion  of  the  Hugoniot 
for  2024-T4  aluminum  indicates  that  the  stress  difference  AC  at  the  inter¬ 
face  amounts  to  about  3  kbar  at  the  stress  level  under  consideration.  The 
offset  DC  at  some  interior  point  in  the  target  would  be  correspondingly 
less.  The  experiments  indicate  that  the  relaxation  of  the  peak  plastic 
wave  stress  from  A  to  E  has  occurred  in  both  hard  and  soft  aluminum  in 
less  than  3.2  mm  of  travel. 

b.  Polyethylene 

The  experimental  data  for  polyethylene  are  presented  in  Table  II. 

Peak  stresses  in  most  of  the  shots  were  about  3  kbar,  but  some  shots  were 
fired  at  about  S  and  11  kbar.  The  depths  of  the  mangamn  gages  from  the 
impact  plane  varied  from  0.78  mm  to  20.3  mm.  No  definite  elastic  precur¬ 
sor  wave  could  be  identified  in  any  of  the  shots.  Pesk  stresses  recorded 
for  any  one  shot  by  the  two  gages  at  different  depths  in  the  target  gen¬ 
erally  agreed  within  the  limits  of  experimental  error.* 

Several  stress  profiles  ar«i  shown  in  Figu  re  9.  A  common  feature  in 
all  the  observed  profiles  was  a  rapid  rise  to  sbout  80  percent  of  the 

*  TV*  ****** i >  g*f«»  *r*  itlurMlb  r«pr*4*ciVI*  u  mkii  1  ptrtnli  Um  ucirttialiti  ii  nWiif  tk* 
p*lr«tkvl*M  r«c«r4*  **r*  f***r«ll>  |r*«t*r  tk**  tk*t. 
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Table  II 


POLYETHYLENE  DATA 


SHOT 

NUMBER 

PROJECTILE 

GAGE 

WAVE 

VELOCITY  DATA* 

Tilt0 

Ui*ec) 

Velocity 

(me//u«ec) 

No. 

Depth 

(mm) 

Streaa^ 

(kbar) 

"o 

(m»/u»ec) 

(mm/^aec) 

ul-2 

(Between 

****»> 

(smAieec) 

B 

D 

12511 

- 

0.256 

- 

19.1 

- 

2.50 

~2.4 

—2. 38 

- 

12512 

0.015 

0.260 

1 

3.22 

1.70 

2.40 

2.04 

2.01 

2 

6.52 

- 

2.24 

2.16 

2.15 

2.29 

12513 

0.005 

0.306 

1 

6.14 

_ 

3.  SO 

2.56 

_ 

2 

12.54 

- 

~2.4 

- 

- 

- 

12514 

0.025 

0.307 

1 

3.115 

2.66 

3.05 

1.98 

1.95 

2 

9.51 

- 

3.00 

2.14 

2.12 

2.23 

12515 

0.037 

0.305 

1 

3.095 

2.65 

3.24 

2.04 

2.01 

2 

6.332 

2.63 

3.13 

2.05 

2.03 

2.05 

1251ft 

0.013 

0.306 

1 

3.194 

2.60 

3. 10 

—2.21 

2.60 

2 

9.634 

2.34 

3.02 

~2.33 

- 

12517 

0.008 

0. 282c 

1 

4.57 

4.29 

5.47 

2.25 

2.22 

2 

10.89 

4.20 

5.41 

2.53 

2.52 

2.79 

12804 

0.022 

0.298 

1 

1.20 

2.  *0 

~3.00 

3.74 

3.50 

2 

20.15 

2.95 

4.20 

2. *>4 

- 

2.98 

12805 

0.022 

0.749 

1 

1.05 

_ 

_ 

_ 

_ 

2 

20.30 

- 

11.0 

2.98 

2.96 

- 

1280ft 

0.011 

0.722 

1 

0.787 

5.0 

9.9 

2.77 

2.58 

2 

18.15 

5.7 

9.8 

2.82 

2.81 

2.82 

* 


Across  6  mm,  the  active  length  of  the  gage. 
See  diagraai  for  definition  of  symbols. 

C  Aluminum  headed  projectile 


3.095mm  THICK 


(a)  SHOT  12,515 


6.332mm  THICK 


4.59mm  THICK 


(b)  SHOT  12,517 


10.89mm  THICK 


0.787mm  THICK  18.15mm  THICK 

(c)  SHOT  12,806 


FIG,  9  MAN  GAN  IN  GAG  6  RECORDS  FOR  POLYETHYLENE 


final  stress,  followed  by  a  much  slower  rise  to  peak  stress.  In  Shot  12,806 
(see  Figure  9c)  one  gage  was  at  a  depth  of  0.78  mm  from  impact  and  the 
other  was  at  a  depth  of  13.1  mm.  Both  gages  recorded  the  same  peak  stress 
but  the  deeper  gage  recorded  a  much  more  rounded  (slower  rising)  pulse. 

The  apparent  double  wave  structure  in  this  record  is  not  being  interpreted 
as  an  elastic  plastic  effect.  The  separation  in  time  of  the  "two  waves" 
is  0.07  fisec  at  0.78  mm  from  impact,  and  only  0.11  ^usec  at  18.1  mm  from 
impact.  The  double  wave  appearance  was  formed  quite  rapidly  and  main¬ 
tained  that  appearance,  with  both  waves  moving  at  nearly  the  same  velocity 
over  a  distance  of  17.3  mm.  The  rounding  off  of  the  top  of  the  wave  front 
profile  with  propagation  distance  has  also  been  observed  in  Plexiglas 
(Ref.  4). 

The  wave  profiles  indicate  a  nonsteady- state  behavior,  and  the  average 
shock  velocities  from  impact  to  the  first  gage  are  lower  in  every  case  but 
one  then  the  average  shock  velocity  between  gages.  Under  nonsteady- state 
conditions  it  is  not  strictly  valid  to  apply  the  Rankine-Hugoniot  jump 
conditions  and  calculate  Hugoniot  points.  However,  assuming  that  the  wave 
has  nearly  reached  steady  state  and  that  the  average  velocity  to  the  final 
gage  is  quite-near  the  instantaneous  velocity,  some  Hugoniot  points  have 
been  calculated.  These  data  are  plotted  in  Figure  10.  Two  points  of 
Wagner,  Waldorf,  and  Louie  (Ref.  5)  are  included  for  comparison. 

In  a  parallel  program  carried  on  simultaneously  with  the  present  work 
[Contract  AF  29(601  ) - 72 14] ,  similar  experiments  have  been  performed  on 
"Teflon."  A  decay  in  peak  atress  of  from  5  percent  to  10  percent  between 
gages  was  observed  (gage  separation  about  6  mm),  as  well  as  a  definite 
rounding  of  the  wave  front  as  it  propagated.  The  wave  profilea  at  a  par¬ 
ticular  gage  location  were  relatively  flat-topped;  thus  the  reduction  of 
peak  pressure  and  change  of  wave  shape  are  not  due  to  norual  attenuation 
but  are  nonsteady- state  effects. 

It  is  of  particular  interest  to  compare  shock  velocities  with  mea¬ 
sured  sound  velocities  in  polymers  because  of  the  dependence  of  sound 
velocity  on  frequency  in  certain  frequency  ranges,  and  because  of  the  lack 
of  evidence  of  an  elastic  precursor  in  any  polymer  yet  studied.  The  data 
for  polyethylene  suggest  that  the  velocity  increased  somewhat  as  the  wave 
propagated  into  the  specimen.  This  effect  could  result  from  a  systematic 
error  in  measuring  the  impact  time,  although  the  effect  was  large  in 
Shot  No.  12. 517  and  a  large  timing  error  would  be  necessary  to  account  for 
it.  The  moat  reliable  transit  times  (and  hence  velocities)  are  those 
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measured  between  like  events — in  this  case,  between  two  gage  records. 

The  corresponding  shock  velocity  at  the  lowest  pressure  of  this  study 
( 3  kbar)  is  2:2  J  0.1  mm/^sec  (three  shots).  This  is  to  be  compared 
with  a  measured  (at  1  Mc/sec)  longitudinal  sound  speed  of  2.26  mm/^sec. 
This  situation  is  in  contrast  to  that  for  Plexiglas,  in  which  the  lowest 
measurements  of  shock  velocity  (~3  kbar)  still  exceeded  measured  values 
of  the  longitudinal  sound  velocity. 

X-ray  patterns  were  taken  of  some  samples  of  polyethylene  before  and 
after  impact.  The  pattern?  prior  to  impact  indicate  a  high  degree  of 
crystallinity  (about  90  percent)  with  random  orientation  of  crystallites. 
Patterns  taken  after  impact  show  that  some  of  the  crystallinity  was  lost, 
but  there  was  a  net  orientation  of  the  remaining  crystallites  such  that 
the  polymer  chains  tend  to  lie  in  a  plane  parallel  to  the  shock  front. 

It  should  be  emphasized  that  these  were  not  controlled  recovery  experi¬ 
ments.  The  specimens  were  subjected  to  gross  deformations  which  can  be 
avoided  by  proper  target  design. 

4.  Summary  and  Suggestions  For  Further  Work  on  Stress  Relaxation 

The  experimental  effort  has  been  directed  toward  the  observation  of 
stress  relaxation  phenomena  in  soft  and  hard  2024  aluminum  and  in  highly 
crystalline  polyethylene.  A  nonsteady- st ate  elastic  precursor  which  de¬ 
creased  in  amplitude  with  propagation  path  length  was  clearly  detected  in 
hardened  aluminum.  Nonst eady - st ate  effects  were  also  detected  in  the  soft 
aluminum  and  polyethylene  but  were  not  as  dramatically  evident  as  in  the 
hardened  aluminum.  A  low- amp  1  it ude  elastic  wave  which  did  not  vary  by  a 
large  amount  was  observed  in  the  soft  aluminum.  No  elastic  precursor  wes 
seen  in  the  polyethylene;  however,  the  wave  shape  was  observed  to  vary  as 
the  wave  propagated,  indicating  a  nonsteady  state. 

The  present  work  represents  only  a  beginning  in  the  much  needed  study 
of  nonsteady- 3tate  propagation  of  low- amp  1 i tude  streaa  waves  and  the  mech¬ 
anisms  of  failure  of  solids  under  rapidly  applied  loads.  An  extension  of 
this  work  could  profitably  be  directed  toward  studies  on  metals  or  poly¬ 
mers.  “Stress  relaxation"  in  the  context  of  this  project  is  manifested  by 
the  decay  in  amplitude  of  the  e'astic  precursor  in  an  elastopiastic  solid 
as  the  wave  propagates  during  the  early  states  of  shock  formation,  t.e., 
very  soon  after  impact.  This  decay  in  aaiplitude  it  associated  with  the 
finite  time  required  for  the  plastic  yielding  to  occur.  The  ultimate 
steady-state  amplitude  of  the  elastic  wave  is  related  to  the  static  yield 
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point  of  the  material.  Inasmuch  as  no  elastic  precursor  has  been  observed 
in  polyethylene  (or  several  other  polymers),  it  is  suggested  that  the  defi¬ 
nition  of  the  term  "stress  relaxation ''  be  broadened  to  include  the  general 
process  by  which  the  wave  approaches  its  steady-state  nrr.file. 

..  a.  Suggestions  for  Work  on  f'otals 

The  present  study  has  established  the  existence  of  stress  relaxation 
in  2024  aluminum.  The  decay  in  amplitude  of  the  elastic  precursor  with 
propagation  distance  appears  greatest  in  hardened  2024  aluminum.  Nonsteady- 
state  behavior  is  visible  up  to  a  depth  of  at  least  0.5  inch  from  the  plane 
of  impact.  It  is  now  clear  that  nonsteady- st ate  effects  are  quite  impor¬ 
tant,  perhaps  equally  as  important  as  nonhydrodynamic  effects,  and  further 
investigation  of  these  effects  is  required  for  a  complete  understanding  of 
wave  propagation  in  solids. 

The  present  work  has  been  devoted  to  observing  the  shape  of  the  pres¬ 
sure  pulse  rise  as  a  function  of  propagation  distance.  The  duration  of  the 
pulses  has  been  essentially  infinite  because  very  thick  projectile  heads 
have  been  used.  Projectile  heads  have  been  made  of  the  same  material  as 
the  target  to  maintain  a  constant  particle  velocity  boundary  condition  an 
the  impact  plane.  Experiments  using  thin-headed  projectiles  should  also  be 
performed,  as  the  long-range  goal  of  this  work  is  to  gain  a  complete  under¬ 
standing  of  shock  propagation.  Such  understanding  must  include  the 
nonsteady- st ate  behavior  of  the  rise  of  the  wave,  the  flow  behind  the  wave, 
and  particularly  the  manner  in  which  rarefactions  overtake  the  shock  and 
produce  attenuation.  At  the  present  time  there  is  every  reason  to  believe 
stress  relaxation  should  affect  all  parts  of  the  flow.  Or.e  "er>  important 
question  to  answer  is  whether  stress  relaxation  affects  th?  late  stages  of 
shock  propagation,  t.e..  wave  shape  after  a  propagation  distance  of  many 
initial  pulse  thicknesaea.  Alto,  what  effect  does  stress  relaxation  have 
on  spall  thresholds  snd  other  damage  criteria'’ 

The  work  suggested  below  is  aimed  at  continuing  the  efforts  of  the 
past  year  toward  answering  these  questions. 

(1)  Aluminum.  Having  established  stress  relaxation  in  aluminum 
using  very  long  pulses,  the  work  must  now  be  extended  to  pulses  of  short 
duration.  Shock  profiles  should  be  observed  after  propagating  several 
pulse  widths  and  results  should  be  compnred  to  pred.ction  of  attenuation 
modela.  The  effect  of  initial  temperature  on  atreaa  relaxation  effects  is 
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unknown.  Experiments  similar  to  those  performed  during  the  past  year 
should  be  performed  on  hardened  aluminum  at  initial  temperatures  from 
200°C  up  toward  melting  (~660°C). 

(2)  Other  Metals.  Very  little  el astoplastir  data  or  stress 
relaxation  data  exist  for  metals  other  than  aluminum  and  iron.  Two 
metals  of  particular  interest  to  which  the  present  work  should  be  ex¬ 
tended  are  tungsten  and  titanium.  Quartz  or  manganin  wire  gages,  as 
appropriate,  should  be  used  to  observe  wave  shape  for  various  thicknesses. 
(The  manganin  wire  pressure  gage  is  particularly  well  suited  for  stress 
propagation  studies  in  polymers.  It  has  already  been  used  quite  success¬ 
fully  in  Teflon  and  polyethylene.)  Initial  conditions  such  as  hardness 
and  method  of  fabrication  should  also  be  varied  as  in  the  case  of  hard 
and  soft  aluminum.  After  investigating  the  behavior  of  long  duration 
pulses,  studies  should  be  made  of  attenuation  with  emphasis  on  late  stage 
behavior  and  spall  prediction  criteria. 

Extension  of  the  present  work  should  first  be  concentrated  on 
further  studies  of  polyethylene.  Such  work  should  include  the  following: 

(a)  Experiments  using  manganin  wire  gages,  as  in  the 
present  work,  should  be  performed  at  higher 
stress  levels  to  determine  the  dependence  upon 
stress  of  the  nonsteady-state  behavior  which  is 
observed  at  low  stresses.  These  experiments 
should  be  perforated  on  both  the  existing  samples 
(for  comparison  with  existing  data)  and  on  newly 
constructed  samples  with  different  microscopic 
structure  ( and  hence  different  mechanical 
properties). 

(h)  The  effect  of  wave  rise  time  upon  shock  struc¬ 
ture  should  be  investigated.  Foam  attenuators 
could  be  used  to  alter  the  shape  of  the  input 
wave.  Manganin  wire  gages  could  then  be  used, 
as  in  the  past,  to  look  at  the  wave  shape  tn  the 
polyethylene  sample  resulting  from  different  in¬ 
put  wave  shapes. 

(c)  The  effect  of  initial  temperature  upon  wave 
propagation  ia  very  important  from  both  the 
applied  and  purely  scientific  viewpoints.  Again, 
bccauae  of  their  temperature  insensitivity, 
manganin  gages  are  very  well  suited  for  the  study 
cf  wave  shapes  in  polyethylene  at  elevated 
temperature!. 


(d)  The  mechanical  properties  of  the  polyethylene 
used  in  the  present  work  appear  to  be  affected 
by  shocking.  The  material  is  more  brittle 
after  shocking  than  it  was  initially  and  X-ray 
patterns  indicate  some  orientation  of  the 
crystallites  which  was  not  initially  present. 

It  is  not  clear  at  present  whether  the  observed 
changes  are  due  to  the  shock  or  to  the  subse¬ 
quent  relief  by  lateral  rarefactions.  To  de¬ 
termine  the  effect  of  shock  on  the  crystal 
structure,  specially  designed  recovery  experi¬ 
ments  should  be  performed  in  which  the  strain 
in  the  sample  is  maintained  oue-dimensional  as 
nearly  as  possible. 

In  addition,  exploratory  work  on  other  plastics  should  be 
carried  out  to  investigate  the  possibility  of  nonsteady-state  effects. 
Carbon  phenolic  would  ,'_.e  a  logical  candidate. 
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SECTION  III 


THEORETICAL  PHASE 

by 

S.  V.  Hanagud  and  R.  C.  Alverson 


1 .  Background 

Research  on  shock  propagation  in  solids  has  been  the  principal 
source  of  information  leading  to  an  understanding  of  the  dynamic  behavior 
and  constitutive  relations  of  these  solids  when  subjected  to  very  large 
stress  and  very  large  stress  rates.  In  the  case  of  metals  in  certain 
stress  ranges,  the  wave  form  consists  of  an  elastic  precursor  followed 
by  a  plastic  shock  wave  which  carries  the  material  to  its  final  stress. 
Some  recent  experimental  work,  (Refs.  6-9)  as  well  as  the  experimental 
work  on  hard  and  soft  aluminum  done  in  the  present  research  program, 
shows  that  the  amplitude  of  the  elastic  wave  is  considerably  higher  than 
that  pedicted  from  elastic-plastic  theories  based  on  static  measurements. 
This  phenomenon  can  be  observed  if  the  measurements  are  made  within  the 
first  few  microseconds  of  the  formation  of  the  two-wave  structure.  As 
the  wave  propagates  into  the  solid,  the  amplitude  of  the  elastic  shock 
wave  decreases  and  asymptotically  approaches  the  amplitude  predicted 
from  the  usual  e 1  as t i c - p 1  as t i c  theory.  The  mechanism  for  such  processes 
of  stress  relaxation  (Ref.  10)  depends  on  the  particular  material.  The 
general  mathematical  description  of  wave  propagation  is  most  conveniently 
given  in  terms  of  Eulerian  coordinates.  In  this  coordinate  system  the 
field  quantities  such  as  displacement,  velocity,  density,  or  stresses 
are  expressed  as  functions  of  the  instantaneous  position  of  a  material 
particle  and  of  the  time.  Using  this  coordinate  system  and  making  the 
assumption  of  uniaxial  strain,  to  coincide  with  experimental  work,  the 
governing  set  of  equations  is: 

da 

,  1 

pv  -  =  0  (conservation  of  momentum)  (l) 

ox 

1  .  3t> 

—  (p)  ♦  —  0  (conservation  of  mass)  (2) 

p  ox 
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(energy  balance) 


(3) 


*  =  Q * —pv 

P 


and 


/(<*  ,  <r  ,  6jP  ^ .  ()  =  0  (constitutive  relations) 


(4)) 


The  symbols  used  in  these  equations  are  explained  as  follows: 


t  -  time 


x  -  Eulerian  coordinate 

B  B 

(•)  =  material  derivative  given  by  ~  ♦  v  -~ 

at  ox 


p(x,t)  * 
v(x, t)  = 
crt(x,t)  = 


density 

particle  velocity 

normal  stress  in  x  direction  (positive  in  tension) 
internal  energy  per  unit  mass 


Q  *  heat  per  unit  mass  added  by  conduction,  radiative 
absorption,  etc. 

£(  =  strain  in  x-direction 

/  *  functions  defining  the  constitutive  relation, 

including  thermodynamic  variables. 


The  formulation  of  appropriate  functions  /  based  on  the  principles 
of  continuum  mechanics  and  the  physical  mechanism  of  the  material  is  the 
first  task.  Later  the  solution  of  the  set  of  partial  differencial  equa¬ 
tions  with  appropriate  initial  and  boundary  conditions  can  be  used  to 
correlate  and  interpret  the  experimental  data. 

In  some  materials,  such  as  polyethylene,  the  effect  of  thermodynamics 
and  viscous  stresses  may  be  more  important  than  the  stress  rate  effects. 
Then  the  functions  /  must  be  defined  to  take  these  into  account. 
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2.  Constitutive  Relations  for  Stress  Relaxation 


The  constitutive  relations  for  stress-relaxing  materials  are  dis¬ 
cussed  here  under  the  assumption  of  uniaxial  strain.  Under  this  assump¬ 
tion  all  dependent  variables  are  functions  of  one  Eulerian  coordinate, 

*,  and  the  time.  Since  there  is  only  one  nonzero  velocity  component, 
the  strain  rates  are  given  by 


3* 


(5) 


and 


•  y 


=  e  =  e 
y  » 


=  0 


The  only  nonvanishing  components  of  the  stress  tensor  are  a..  v 
and  cr t,  and  because  of  symmetry  under  uniaxial  3train  we  have  a,  *  crj. 
Decomposition  of  the  stress  tensors  into  a  hydrostatic  pressure  p  and 
stress  deviator  f)  gives 
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(7) 


The  formulation  of  the  constitutive  relation  consists  of  determining 
the  relation  between  p,  DtJ>  «nd  the  kinematic  variables.  It  will  be 
assumed  that  the  hydrostatic  pressure  p  is  a  function  of  density 
alone,  t.e., 

P  ’  fip)  (8) 
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The  function  f(p )  depends  on  the  particular  material.  Two  suitable  forms 
of  f(p)  are 


When  the  function  f(p)  is  specified  the  remaining  task  is  to  specify 
the  relationship  between  the  stress  deviator  and  the  kinematic  variables. 

a.  Elastic-Plastic  Model 

The  very  commonly  used  relation  between  the  stress  deviator  and  the 
kinematic  variables  is  the  familiar  elastic-plastic  relation.  This  can 
be  described  as  follows.  In  the  case  of  uniaxial  strain  there  is  only 
one  iruependent  component  of  the  stress  deviator  because 


'  -2°,,  '  -20..  '  JK  -  V  •  <»> 

In  a  similar  wi)  the  components  of  the  strain  deviator  Et)  are 


£  -  --He  >  (10) 

. .  3  * 


The  strain-rate  deviatora  ire 
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If  the  solid  is  defo-juning  elastically,  i.e.,  if 


Id, J  1  y.  y  lo 


then 


2GE. 


(ID 


(12) 


(13) 


In  Eqs.  (12)  and  (131,  y  is  the  deviatoric  yield  stress*  and  G  is  the 
modulus  of  rigidity.  Furthermore,  Tresca  type  of  yielding  has  been 
assumed.  The  solid  is  said  to  be  deforming  plastically  in  regions  where 

ID.J  5  y  ■  (14) 

Once  |Dif|  has  attained  the  yielding  value  of  y,  | |  cannot  increase 
beyond  this  value.  However,  the  strain  deviator  can  continue  to  increase. 


b.  F.last  ic-Plastic  Relaxing  Model 

The  elast ic -plast ic  model  cannot  account  for  certain  experimental 
observations  in  shock  wave  propagation.  Measurements  made  within  the 
first  few  microseconds  have  shown  that  the  two-wave  structure  predicted 
by  the  e  1  aa t ic -plat t ic  model  ia  not  formed  instantaneously  but  is  only 
realised  after  some  time  by  a  process  of  relaxation.  During  the  relax¬ 
ation  process,  the  amplitude  of  the  elastic  wave  (or  elastic  precursor) 
ia  considerably  higher  than  that  predicted  from  the  elastic-plastic 
model.  As  the  wave  propagates  into  the  solid  the  amplitude  of  the  wave 
decreases  and  asymptotically  approaches  a  value  corresponding  to  the 
static  yield  stress. 

*  J  *  2 ' 3 f  tfctrt  T  It  ik»  yit|4  ttrttt  i*  tioplt  itttitt. 


33 


(1)  Ma  croscopic  Relaxation  Equations 


The  elastic'plastic  stress-strain  relationships  outlined  above 
are  macroscopic.  They  can  be  generalized  to  include  stress  relaxation. 

A  simple  model  can  be  constructed  by  considering  the  rate  type  of  the 
constitutive  relation  in  the  form 

D  »  F(D,E,p )  (15) 

where  E  is  the  rate  of  deformation  tensor.  The  symbol  A  represents  the 
Janmann  or  co-rotational  type  of  "tress  rate  (Ref.  11).  In  the  case  of 
uniaxial  strain  these  stress  rates  reduce  to  the  usual  time  derivatives 
along  particle  paths. 

A  speci fic  relaxing  type  of  stress-strain  relationship  as  applied 
to  the  case  of  uniaxial  strain  can  be  written  as 


if 


otherwise 


£ 


I  M 


2  GT 


>  y 


(16) 


(17) 


(18) 


In  these  equations  8  can  take  the  value  of  i  1.  This  accounts  for  the 
fact  that  yielding  can  take  place  when  the  absolute  value  of  the  deviator 
reaches  y.  T  can  be  called  the  relaxation  constant  and  has  the  dimension 
of  time. 


The  first  term  on  the  righthand  side  of  Eq.  (16)  takes  care  of 
the  fact  that  the  components  of  the  stress  deviator  can  hsve#vaiues 
greater  than  the  static  deviatoi ic  yield  stress.  This  then  can  result  in 
the  amplitude  of  the  elastic  precursor  wave  being  much  larger  than  the 
equilibrium  value.  The  form  of  the  term  suggests  thst  the  devistor 
relaxes  with  time  towards  the  static  yield  value.  The  second  term  on  the 
righthand  side  of  Eq.  (16)  takes  care  of  the  elastic  deformation.  This 
equation  can  be  written  as 
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(19) 


“  -f  *  2  GE 
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with  the  condition  that 


D  =  2G£ 

*  *  *« 


(20) 


when 


Id 


»  y 

XX  J 


(21) 


Alternatively  it*can  be  stated  that  when  EtI  *  0  and  D  *  8y\  then 


D  -  0 

X  X 


(22) 


If  Eq.  (19)  is  expressed  as 


D  ,  •  F(D  ,E  ,p) 

XX  XJ1  u,r 


(23) 


we  have 


In  a  similar  way 


F(0y,O,p)  -  0 


D  -  F(D  ,E  ,p) 

yy  yy'  yy'^ 


*  no„.E,ttp) 


(24) 

(25) 


However,  these  equations  will  be  identical  to  Eq.  (19)  in  case. of  uniaxial 
strain;  y  may  be  a  function  of  some  invariants  of  the  stress  tensor  in  the 
general  tensor. al  constitutive  equation. 

Further  generalisations  can  be  made  by  considering  general  forma 
of  functions  F  in  Eq.  (19).  However,  these  generalisations  should  be 
based  on  the  physical  mechanisms  and  experimental  observations. 

In  the  case  of  hard  and  soft  aluminum,  experimental  observation 
indicates  that  a  two-wave  structure  is  formed  a  very  short  time  after 
impact.  However,  the  amplitude  of  the  elastic  precursor  is  initially  very 
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large  compared  with  the  value  predicted  by  the  elastic-plastic  model. 

The  elastic  wave  amplitude  decreases  within  a  few  microseconds  to  the 
steady-state  amplitude  of  the  elastic  precursor  as  predicted  by  the 
elastic-plastic  theory.  This  behavior  indicates  that  the  yield  stress 
is  very  large  at  early  times,  relaxing  to  the  static  value  in  a  few 
microseconds.  (This  can  be  called  the  yield  relaxation  time.)  It  also 
indicates  that  the  relaxation  of  the  elastic  precursor  towards  the  equi¬ 
librium  value  is  taking  place  much  faster  than  the  yield  relaxation  time 
and  thus  is  setting  up  the  two-wave  structure  at  a  very  early  time.  A 
model  to  explain  this  behavior  can  be  obtained  by  adding  to  Eq.  (16)  an 
additional  equation, 


y  >  y  o 


In  this  case  y  is  no  longer  a  constant  but 

y  =  y(crt,cry,p) 

and  decreases  with  time  towards  the  static-yield  value 

y  5  >0 


(26) 


Another  way  of  generalizing  the  relaxation  equation  is  as 
follows.  Equation  (16)  and  Eq .  (8)  can  be  simplif'ed  by  substituting  in 
them  the  expressions  fo;  p  and  Dtl  in  terms  of  o ^ ,  and  eliminating  c 
between  the  two  equations.  Then  we  have 


where 


end 


♦  F(p)p  ■ 


F{p> 


AC,  df_ 

3  p  dp 


F'(p,a) 


fip)  G‘  9y0 
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(27) 
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Alternatively  one  can  think  of  Eq.  (27)  as  a  one-dimensional  relaxation 
equation  with 


F(p) 


da 

f 

dp 


the  quantity  c2  is  the  square  of  the  sound  speed  and  F'(p,a )  is  the  dis¬ 
sipation  function.  Then  in  any  one-dimensional  problem  a  function 
a x  •  ~F(p)  and  F'(p,a*)  should  be  specified.  a*  =  ~F(p)  will  be  the 

elastic  a*t  p  relationship.  This  can  be  found  from  many  steady-state 
experimental  results.  A  function  for  F'(p,a )  based  on  Zener's  model  is 


F' (p,at) 


P. 


where  p (  is  the  amplitude  of  elastic  precursor  when  complete  relaxation 
has  taken  place.  This  would  be  the  amplitude  of  the  elastic  precursor 
on  the  basis  of  elastic-plastic  models.  Another  form  for  F' (p,a )  can 
be 


F'ip.a) 


P,  -  P,(CT,.P) 


and 


P,  '  P 

P,  - - = -  (28) 

'  2 

where  Tj  is  the  relaxation  constant  for  the  yield  stress. 

The  value  of  p(  on  any  particle  is  of  significance  from  the 
time  the  particle  comes  under  the  influence  of  the  leading  elastic  ahock 
wave.  If  this  time  u  specified  as  t(  and  the  pressure  at  thia  instance 
is  p  ,  then 

P,(*.>  ’  7<P,> 

This  equation  specifies  the  iuitial  condition  for  Eq.  (28).  A  similar 
reasoning  can  be  used  to  specify  the  initial  condition  for  Eq.  (26). 
Calculation  of  wave  profile  with  these  equations  will  be  discussed  later. 
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(2)  Micrcacopic  Relaxation  Equations  (Dislocation  Mechanism) 


Shear  yielding  and  resultant  plastic  flow  in  metals  are  due  to 
the  motion  of  dislocations.  A  sufficient  number  of  mobile  dislocations 
is  necessary  to  cause  yielding.  In  dynamic  problems  a  finite  tirrnt  ,.s 
necessary  for  dislocation  production  and  movement  and  hence  for  production 
of  plastic  flow.  This  finite  time  is  then  the  cause  of  the  large  initial 
amplitudes  of  elastic  precursor  waves. 

It  is  known  that  a  sufficiently  strong  external  stress  can 
cause  a  dislocation  to  unpin  itself  and  hence  become  mobile.  Dislocations 
can  also  be  created  by  the  Frank-Read  mechanism,  a  kind  of  multiplication 
or  breeding  process,  and  again  the  effective  agency  is  stress.  Once 
unpinned,  a  dislocation  is  free  to  move  and  contribute  to  yielding  or 
plastic  strain,  until  it  again  becomes  trapped  at  another  grain  boundary 
or  impurity.  The  total  number  of  moving  dislocations  effectively  con¬ 
tributing  to  the  plastic  strain  rate  «t  any  moment  is  the  net  result  of 
competition  between  the  processes  of  unpinning  or  creation  (sources)  and 
of  trapping  (sinks). 

It  is  natural  to  suppose  that  unpinning  is  an  excitation 
process  with  a  rate  depending  on  the  excess  of  shear  stress  over  yield 
stress,  while  the  Frank-Read  mechanism  produces  new  dislocations  from 
existing  ones  at  a  rate  that  is  proportional  to  the  number  of  existing 
ones  and  that  also  depends  on  the  excess  shear  stress.  To  simplify  the 
analysis  we  approximate  the  excitation  process  by  taking  the  rate  simply 
proportional  to  the  excess  shear  stress.  The  rate  of  crapping  is  of 
course  proportional  to  the  number  of  dislocations  present.  These  rela¬ 
tionships  can  be  put  in  mathematical  form  in  the  following  way.  Suppose 
that  at  some  particular  time  the  mobile  dislocation  density  is  Nm  «nd  the 
total  dislocation  density  is  N(,  hence,  the  density  of  pinned  dislocations 
is  (A,  -  Nm ) .  Assume  that  the  pinned  dislocations  become  Frank-Read 
sources,  ao  that  a  term  g(Z(A(  -  Nm )  contributes  to  dN^dt.  The  constant 
g,  represents  the  strength  of  the  Frank-Read  source  and  Z  is  defined  to 
be  (t  -  E/2>,  the  excess  shear  stress  above  the  yield  point.*  The  pinned 
dislocations  may  also  become  mobile  by  various  freeing  mechanisms  and 
this  contributes  another  term,  g2 Z{Nf  -  Nm ) ,  to  the  rate  dNm/dt.  ie 
assume  further  that  mobile  dislocations  can  become  pinned  by  interaction 


f  ■  1/2  y  .  ■'  ’  1/lUr  -  e  >;  k*»c*  t  •  J/4tf  - 

»  y  »*  7 


31 


with  impurities,  etc.,  thus  tontributing  a  rate  term  of  the  fo.m 
-( 1 /Ty )Nm (v/vm ) ,  where  v  is  the  dislocation  velocity  and  its  maximum 

value.  Mobile  dislocations  can  also  simply  disappear  irreversibly  by 
migration  to  free  surfaces  and  voids,  giving  a  similar  t erm - ( 1 / T  2 )Nm ( v/vm ) . 
These  assumptions  result  in  the  following  expressions  for  the  rate  of 
change  of  the  mobile  and  total  dislocation  densities: 


dN 


4  /II 

»  — g Z(N,  -  N  )  -  (  —  ♦  —  \  N  I  — 
dt  3  8  '  -  ‘  T  r  1  -  ! 


(r.  *  0*-  (”■) 


(29) 


Here  g  »  (3/4)(gj  ♦  g 2 ) ,  and  e  *  gl/(gl  *  g2).  It  can  be  seen  from  the 
definition  that  e  is  a  number  between  vero  and  unity.  When  £  is  close  to 
rero,  a  pinned  dislocation  is  more  likely  to  become  freed  than  to  create 
a  new  dislocation.  A  value  of  £  close  to  unity  means  that  a  pinned  dis- 
location  will  probably  create  a  new  one  rather  tha"  become  freed.  If  the 
stress  history  at  some  point  in  the  solid  is  known,  Eq.  (29)  can  be 
integrated  to  give  the  total  and  mobile  dislocation  densities  at  this 
point. 

In  order  to  get  some  idea  of  the  choice  of  parameters  in  this 
model,  we  consider  the  application  of  a  constant  stress.  The  ratio 
(r  ■' i>  )  will  then  be  constant,  since  the  dislocation  velocity  is  assumed 
to  be  a  function  of  the  -tress  alone.  Eliminating  from  Eq.  (29)  gives 


d'S 


dt 


(  I 


(0(r,  ■  0 


dt 


*  ?•*  r 


(i  - 1) 


t 

F 


We  want  the  mobile  dislocation  denaity  to  irrreese  with  the  appl. cation 
of  stress.  To  ensure  this,  the  equation  above  must  have  a  solution 
which  increases  with  time.  This  requires  that 
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(1-6)  £ 

-  -  —  <  o 

T2  Tl 

This  condition  will  be  satisfied  if  we  choose  e  close  to  unity  and 

T2  >  7V 

The  dislocation  velocity  is  assumed  to  depend  on  the  applied 
shear  stress.  Following  Gilman  and  Johnston-  (Ref.  12). 


Now  the  shear  strain  r t  is  given  by 

y  =  bvS 

'  ■ 

where  b  is  the  Burger's  vector.  The  remaining  task  is  to  interpret  y 
in  tcms  of  e  . 

X 

Under  the  conditions  of  unicxial  strain 


€  *6-0 

r  « 

Furthermore,  the  total  strain  can  be  split  i  ito  elastic  and 
plastic  strains.  Then 
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Because  of  symmetry 


Furthermore 


s'  ♦ 

X  X 


■  e'  *  s'" 
y  y 


c '  +  e? 

I  2 


e*  +  eP 
y  y 


If  it  is  assumed  that  the  plastic  strain  does  not  contribute 
to  any  volume  change 


+  2ei>  =  0 

*  y 


Under  uniaxial  strain  the  maximum  shear  strain  is  oriented  at  45°  to  the 
x  axis  and  the  maximum  plastic  shear  strain  rate  is  given  by 

1  . 

y  =  r(e;  -  eM 


The  relaxation  function  is 


F  -  2Ge  p 


8 

TGbA 


where  6  is  the  Burger’s  vector  and  Nm  must  be  calculated  from  Eq.  (29). 
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c .  Thermodynamic  Considerations 

The  results  of  many  wave  calculations  are  in  disagreement  with 
experiment.  There  are  many  reasons  for  these  discrepancies,  one  of 
which  is  undoubtedly  the  fact  that  the  calculations  r.re  based  on  consti- 
tt  ive  relations  which  are  independent  of  tit..e;  i.e.,  stress  relaxation 
and  viscous  effects  are  neglected.  However,  in  some  materials  the  mac¬ 
roscopic  properties  may  be  highly  temperature-dependent,  and  theoretical 
predictions  should  therefore  be  based  on  thermodynamic  as  well  as  purely 
mechanical  considerations.  Most  theories  of  wave  propagation  in  solids 
are  based  on  a  governing  set  of  differential  equations  which  is  indepen¬ 
dent  of  energy  considerations.  It  is  the  purpose  in  this  section  to 
show  how  a  theory  which  includes  thermodynamic  considerations  may  be 
constructed  for  a  stress-relaxing  solid.  We  consider  the  perfectly 
elastic,  perfectly  plastic,  and  the  plastic  stress- relaxing  solids,  in 
that  order. 

If  we  let  T,  S,  and  E  denote  the  temperature,  entropy,  and  internal 
energy  per  unit  mass,  respectively,  we  have  from  the  first  and  second 
laws  of  thermodynamics 


TdS  *  dE  -  a  ijde 


where  a .  are  the 
and  we  assume  the 

a  e  s  a  e 

i i  . i  ii  ii 

in  the  form 


stress  components  and  are  the  strain  components 

summation  convention  over  the  repeated  subscripts: 

*  a,y  €xy  *  etc.  The  relation  above  can  be  rewritten 


» J 


.'he  constitutive  relations  considered  in  the  present  work  for 
purely  elastic  behavior  in  the  case  of  uniaxial  strain  are 


and 
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These  equations  can  be  solved  for  cr  and  a  in  turn  and  the  results 
can  be  substituted  into  the  thermodynamic  relation  ebove.  It  is  thus 
readily  seen  that  if  a  t  and  a are  independent  of  T,  that  is,  if  G  and 
f(p)  are  independent  of  temperature,  then 


and 


do 

X 

17 
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dE 

de 
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o .  , 
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This  equation  can  be  integrated  directly  to  give  the  internal  energy, 
once  the  strains  are  given.  Thus  the  purely  elastic  model  which  is 
mechanically  consistent  does  not  violate  either  the  first  or  the  second 
law  of  thermodynamics.  If,  in  the  constitutive  relations  given  above, 
either  or  both  G  and  /  depend  upon  the  temperature,  we  then  append  to 
this  system  of  equations  the  equation  which  expresses  conservation  of 
ene  rgy : 


dQ  dE  _  1 

17  =  p  7t  1 ai> 


In  this  equation  Q  is  the  heat  flux  or  heat  input.  This  equation, 
together  with  the  constitutive  relations  and  the  equations  governing 
conservation  of  mass  and  momentum  and  the  first  and  second  laws  of 
thermodynamics  then  give  a  consistent  set  of  equations  with  which  theo¬ 
retical  studies  of  wave  propagation  in  a  perfectly  elastic  heat-conducting 
medium  can  be  made. 

We  now  consider  the  perfectly  plastic  solid  for  simplicity,  and  use 
the  von  Mises  yield  criterion 


e 
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where  e  is  the  strain  rate  and  D  is  the  stress  devictor. 

ij  i ) 

manipulation  of  the  above  we  have 
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We  also  assume  here  for  simplicity  that  the  material  is  incompressible. 

The  equation  of  heat  conduction  or  conservation  of  energy  then  becomes 

dQ  dE 

Tt  ‘  It 

Again,  if  the  yield  stress  is  independent  of  temperature  the  mechan¬ 
ical  problem  is  separate  from  the  thermal  problem  and  the  latter  may  be 
solved  after  the  mechanical  problem  has  been  solved. 

The  case  of  a  stress  relating  elastic-plastic  medium  is  much  more 
complex  but  a  combination  of  the  procedures  given  in  the  two  cases  above 
applies.  The  system  of  equations  which  is  both  thermodynamically  and 
mechanically  consistent  is  given  by  the  conservation  laws  'mass,  momentum, 
and  energy),  the  constitutive  relations,  and  the  combined  expressions  for 
the  first  and  second  laws  of  thermodynamics.  If  the  material  properties 
are  independent  of  the  temperature  the  problems  are  separate,  but  if  f(p ) 
and/or  the  relaxation  constants  depend  on  the  temperature,  then  the  heat 
conduction  problem  must  be  solved  simultaneously  with  the  mechanical 
problem. 

These  considerations  have  been  based  on  the  ideas  of  equilibrium 
thermodynamics,  and  to  date  no  computations  have  berfn  attempted  using  the 
theory  discussed  above.  For  simple  cases  such  computations  are  entirely 
feasible  if  sufficient  data  on  the  temperature  dependence  of  the  consti¬ 
tutive  relations  are  known.  Further  extension  of  this  work  to  nonequi¬ 
librium  thermodynamics  depends  on  formulating  a  definition  of  entropy 
which  is  consistent  with  the  equilibrium  case.  The  thermodynamic  con¬ 
siderations  are  likely  to  be  important  in  some  materials,  since  high 
temperatures  may  be  obtained  at  impact  surfaces  and  the  resulting  tem¬ 
perature  distribution  may  then  alter  the  wave  shape  and  may  be  a  source 
of  high  energy  dissipation. 
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3.  Method  of  Computation 


a.  Characteristics  of  the  Differential  Equations 

The  equations  governing  the  wave  propagation  under  uniaxial  strain 
can  be  written  in  the  following  way. 

Equation  expressing  conservation  of  mass: 


Bp  3  p  B  v 

_  +  V  ♦  n  — 

3 1  3 x 


Newton ’ s  i aw: 
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(31) 


(32) 


In  these  equations  x  is  the  Eulerian  coordinate,  a  *  is  the  stress  in  the 
jc-direction,  v  is  the  particle  velocity  in  the  x-direction,  an i  P  is  the 
density.  Equations  (31)  ' nd  (32),  together  with  the  appropriate  consti¬ 
tutive  relationships,  complete  the  governing  equations. 

Constitutive  relations: 


(33) 
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In  these  equations  the  dots  represent  the  derivative  along  the  particle 
paths.  Furthermore,  Eq.  (34)  assumes  the  simple  one- pa r ameter  macroscopic 
model  for  the  stress  relaxation.  This  equation  will  be  simplified  and  put 
in  a  form  such  that  it  can  be  generalized  to  include  other  moJels. 

By  using  the  following  relationship 


Equation  (34)  becomes 
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where 


>  "  2  y 


By  definition 
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Then  from  Eq.  (31) 
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Equation  (35)  now  becomes 
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This  equation  and  Eq.  (33)  can  be  used  to  eliminate  cr^.  Then  the  follow¬ 
ing  equation  is  obtained; 


-|4£  ♦  ii)b  f(P) 
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or 


wher; 
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(41) 
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T  T  T  3 


(42) 


Equation  (40)  is  in  a  form  in  which  most  of  the  constitutive  relationships 
involving  stress  relaxation  and  only  the  first-time  derivative  of  stress 
can  be  written.  Equation  (40)  can  also  be  written  as 


Then 


da  t  ■  -F(p)dp  -  F'  \.p,a  t)dt 

dp  F'(p,aM)  dat  j 

dt  "  '  f(p)  '  dt  F( p) 


(43) 

(44) 
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The  equation  expressing  the  conservation  of  mass  can  now  be  rewritten  as 


“a»  ^ 

— —  ♦  F'ip.cr  )  •  pF{p)  — 

C It  cx 


By  introducing  the  notation 


P*  ' 


equation  (45)  becomes 


dP,  dv 

-  ♦  F'(p,cr  )  ”  pF(p)  - 


The  use  of  the  definition  of  the  derivative  along  the  particle  path  yields 


3P,  3P, 


Pc'1  -r-  =  F'(p,pt) 


The  equa  ion  oi  motion  [E<«.  (32)]  can  also  be  written  as 


dv  dt; 


c  u  +  tx  +  p  irf  =  0  • 


By  adding  Eqs.  (48)  and  (49)  we  have 


v  dP.  ['dv 

— —  ♦  (u  *  c)  *  pc ■—  *  (v  ♦  c) 

dt  dx  r  dt 


Similarly,  by  subtracting  Eq.  (49)  from  Eq.  (48) 


.  dP.  fdv  dr 

♦  (v  -  c)  — - Pci—  ♦  (r  -  c)  -r- 

dt  dx  dt  d* 


The  form  of  Eqs.  (50)  and  (51)  suggests  that  they  can  be  written  as 


"3T  • 
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and 


where 


(53) 


represent  the  total  derivatives  along  the  paths  Ot  and  yS  defi,  ed  by 


and 


dx 

dt 


v  *  c 


(54) 


/?: 


dx 

—  -  V  ~  c 

dt 


UJ) 


The  paths  a  and  /?  define  the  chara.  eristics  (Ref.  13)  of  the  set  of 
equations  governing  the  probler  ar  .  Eqs.  (52)  and  (53  1  are  the  two  com¬ 
patibility  relations  valid  along  ach  of  these  cha rac t <  r is ' ic s.  Thes. 
equation  •  can  be  used  to  set  up  a  sc  erne  of  numerical  integration  oi  t.l  » 
initial-boundary  value  problem. 

These  equations  are  valid  only  if  |0  |  >  v .  If  j  D  f  ^  j  £  y,  the 

material  is  in  the  elastic  regime.  Equations  (54)  and  (55)  for  the  har- 
acteristics  are  still  valid;  however,  the  compatibility  relations  are 


An  alternate  way  of  deriving  the  equations  of  characteristics 
and  the  compatibility  relations  based  on  the  existence  proof  of  Csuch  - 
Kowalewski  is  explained  in  Appendix  II. 
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b.  Numerical  Method 


The  mesh  points  are  taken  to  be  fixed.  The  fixed  interval  Ax  is 
along  the  x-axis,  and  At  along  the  time  axis.  The  ratio  Ax/At  is  deter¬ 
mined  from  a  knowledge  of  the  maximum  values  of  particle  velocity  and 
the  sound  speed  c.  The  set  of  fixed  mesh  points  is  shown  in  Figure  11. 
The  curve  labeled  I  represents  the  interior  boundary  and  the  curve 
labeled  II  represents  the  leading  shock  wave  penetrating  into  the  undis¬ 
turbed  medium. 


FIG.  11  FINITE  DIFFERENCE  MESH  FOR  THE  METHOD  OF  CHARACTERISTICS 

The  numerical  determination  of  the  field  quantities  at  a  mesh  point 
x  =  x0+nAx,  t  =  jAt  differs  according  to  whether  the  point  is  on  either 

Curve  I  or  II  or  on  an  interior  point,  in  Figure  11. 

(1)  Numerical  Determination  at  a  General  Interior  Point 

In  Figure  12,  it  is  assumed  that  the  point  labeled  x  lies 
somewhere  on  a  horizontal  line  t  -  t  j,  and  that  the  points  Xj,  x2,  and  x3 
are  on  a  horizontal  line  t  ~  t j  “  A{ ;  the  x  coordinates  of  x2  and  x3 
are  x,  -  Ax  and  Xj  +  Ax,  respectively. 

It  is  also  assumed  that  the  values  of  pj ,  p,  and  v  are  known  at 
points  Xj,  x j ,  and  x3<  x2  must  lie  either  on  or  to  the  right  of  Curve  I 
oi  figure  11  and  x3  must  lie  either  on  or  to  the  left  of  Curve  II.  The 

values  of  the  field  quantities  at  x  are  to  be  determined. 
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FIG.  12  TYPICAL  MESH  POINT  FOR 
A  GENERAL  INTERIOR  POINT 


The  f i  rst  step  is  to  ap¬ 
proximate  the  characteristic  curves 
through  x  by  straight  lines  an  1  find 
their  intersections  with  the  hori- 
lontal  line  t  »  tj  -  At.  To  this  end 
we  approximate  the  curve  with  slope 
dx/dt  *  v*  by  the  straight  line 
through  x  with  slope 


dx 

dt 


i  (  x  , )  ♦  v  (  x  2  ) 


2 


(56) 


This  curve  represents  the  particle  path  along  which  fcq .  (40)  is  valid. 

The  intersection  of  this  line  with  t  ■  tj  -  At  is  then  denoted 
by  in  Figure  12.  The  x  coordinate  of  x|  is  then  readily  computed 
through  the  formula 

tv( Xj)  ♦  v( Xj)  ] 

*1  - - T - At  ♦  *,  •  (57) 


From  Eq.  (57)  we  note  that  x'j  may  be  greater  than,  equal  to, 
or  less  than  x,,  according  to  whether 


v( x  t )  ♦  t( x,  )  > 

■  9 

2  < 


0 


If  x  |  X|  we  de  f  i  ne 


and 


if 


(58) 


(59) 


Tb«  i«m  Wi  4t  *  *  rarimli  tbt  ilurK'fnitu  tint  • 
•  •lid  <  •*•  Appcadu**  It  *>d  III). 


leaf  »bicb  lb* 


r*l*UM»  *r* 


The  values  of  px,  p,  and  v  at  x}  are  then  found  by  interpolation 
through  the  formula 

F<*;>  =  F(A)  ♦  kjtFfB)  -  F(A)] 

where  F  is  successively  equal  to  px,  p,  or  v,  and  A  is  x1#  B  is  x2  if  Xj  is 
given  by  Eq.  (58),  and  4  is  xx  and  B  is  x3  if  is  given  by  Eq.  (59). 

•In  a  similar  manner  points  x2  and  x3  of  Figure  11  are  found  as 
the  intersections  of  the  straight  lines  which  approximate  the  curves 

dx 

-  «=  \)  i  C 

dt 


Thus 


t  *  t  j  -  At 


[( v  *  c)(xj)  *  (v  +  c)(x2)] 


At  ♦  x3 


[(u  “  c)(Xj)  +  iv  -  c)( X3)l 


At  ♦  x. 


In  the  equation  above,  (v  ♦  c)(Xj),  etc.  denote  the  value  of 
(v  ♦  c)  at  point  x3.  It  is  to  be  noted  that  x2  lies  between  x2  and  Xj, 
and  x3  lies  between  x3  and  x3  because  of  the  way  in  which  the  ratio  Ax/A t 
was  taken.  We  then  define 


•2  x  2 


x  3  “  *1 


and  the  values  of  p  x ,  p,  a::J  -  at  x'2  and  x'3  are  then  given  by  Eq.  (60) 

F(x;>  =  F(A)  *  \j [F(B)  -  F(A)] 

where  x[  is  successively  x2  and  x3,  etc. 
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The  finite  difference  equations  which  must  be  solved  to  give  the 


of  pjt 

p ,  and 

v  at  a.-  can  be  written  in  the 

fol 1  owing  way. 

P  Jx)  - 

P  Jx'2) 

♦  /  ( *2 ) ( X 2 )  [f ( * )  ”  v( X 2 ) 1 

=  F'lx'2)JSt 

(65) 

Px(x)  -  pt(x‘3 ) 

-  p(  X  3  )  C  (  X  3  )  [  V  (  Z  )  -  !)  (  t  J  )  ] 

-  F ' ( x 3 )At 

(66) 

In  these  equations 

F‘  (x[)  =  F‘  [Mx’2)  ,pz(x!i)} 

F'ix  3)  =  F'  [plx‘z)  ,p,(x'3)] 

Px(x)  -  pjx[)  =  -F(x\)[p(x)  -  fAx[)\  -  F'(x[)AT  . 

First  p  x(  x)  ,  ti(r),  and  p(  x)  are  computed  from  these  equations  with 
F'  =  0.  Then  D  is  calculated  from  the  following  equation 


4 G  p(x) 

—  log  - 

3 


1  1  2  - 

\D,X\  <Ty 


the  calculations  are  complete.  If 
1  I  2  ~ 

\D  >  ~  y 

Eqs.  (65)  to  (67)  are  used  to  compute  px(x), 
v( x) ,  and  pi  x)  with  F'  in  all  the  equations. 
In  these  calculations  the  value  of  6  is  +1  if 
D  >  0,  and  6  =  -1  if  !)  <  0. 

( 2 )  Numerical  Determination  at  a 
Point  on  the  Boundary 

In  this  case  the  point  x  lies  on 
Curve  I  and  is  not  necessarily  one  of  the 
a  priori  fixed  mesh  points.  Referring  to 
Figure  13,  the  values  of  the  field  quantities 
are  known  at  x)(  x 3,  xt,  and  they  must  be 


^1  •>!  •«)  •t4 

M-fTM-7 

FIG.  13  TYPICAL  POINT 
ON  THE  INTERIOR 
BOUNDARY 
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determined  at  x.  Note  that  * ,  and  x  need  not  be  mesh  points  in  the  fixed 
mesh  but  are  both  points  on  Curve  ’  >d  hence  on  a  particle  path.  The 

point  x  is  then  found  from 

x  c  r(x, )A t  *  x j  171) 

where  it  has  again  been  assumed  that  point  x  lies  on  t  =  t ,  and  points 
Xj,  xj,  x4  lie  on  t  =  f,  "At. 


The  point  x'  is  then  located  by 

[U  -  c)(x.)  *  (v  -  c)(x3)l 


Two  cases  can  arise 


in  tq.  (72): 


and 


If  Eq.  (73)  is  satisfied,  then 

x  3  -  x  , 


and  if  Eq.  (74)  i"  satisfied 


(73) 

(74) 

(75) 

(76) 


In  either  case,  the  values  of  the  field  quantities  at  x3  are 
found  by  making  the  proper  but  obvious  substitutions  in  tq.  (60). 

The  finite  difference  equations  which  must  be  solved  to  yield 
the  field  quantities  are  then: 


(*)  -  P  (  X  3 1  -  p(  x'  )  c(  X  j)  U(  x)  -  f(*3)J 


pjx)  -  P^Xj) 


-F(x,)fp(x)  -  p(xj)]  -  F'(x,)At 


(78) 


Along  Curve  I,  pt(x)  o?  v(x)  is  known.  Then  Eqs.  (77)  and  (78)  can  be 
used  to  determine  the  remaining  two  of  the  three  quantities  pf(x),  v(x), 
p( x) .  Again,  the  quantities  are  computed  with  F'  =  0.  The  stress  devia¬ 
tor  D  is  computed  from  Eq.  (78).  Then  if 


Id  !  <  -  y 

'  *« '  3  7 


the  calculations  are  complete.  Otherwise,  the  calculations  are  done 
using  Eqs.  (77)  and  (78)  with  F' .  Again  6  =  *  1  if  Dti  >  0  and  6  =  -1  if 

D, ,  <  0. 

(3)  Numerical  Determination 
at  a  Point  on  the  Shock 
Wave 


•*4 


We  now  consider  the  point  *  to 
be  on  the  shock  wave  as  shown  in  Figure  14 
In  this  instance  the  points  x  and  x.  lie 
on  the  shock  wave  and  are  not  in  general 
fixed  mesh  points.  The  coordinate  of  x  is 
given  by 


FIG.  14  TYPICAL  POINT  ON 
THE  SHOCK  WAVE 


U(xl ) At  ♦  *, 


(79) 


where  U(x j)  is  the  shock  wave  velocity 
at  Jtj.  The  coordinate  of  x'2  is  given  by 


[ ( v  *  c)(x 2 )  ♦  (v  ♦  c)(*j)] 


At  ♦ 


(80) 


Again  in  Eq.  (80)  two  cases  can  arise:  If  x2  >  x2 


and  if  x'2  <  x2 


Ax 
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The  values  of  the  field  quantities  at  x'2  are  found  by  applica¬ 
tion  of  Eq.  (60'.  The  field  auantities  at  *  are  then  determined  by  the 
equa  t ions 


P0U(x)  = 

p(  x)  [t/(  x)  "  v(  x)  ] 

(81) 

P„<*>  ■ 

Pq  U ( x) v ( x) 

(82) 

p(xj)c(xj)  [  i>  ( x ) 

-  V  (  X  2  )  3 

F' ( x  j )Aj 

(83) 

x  4G 

*)  =  *  V  lo£ 

p(x) 

- —  *  /(p) 

(84) 

Po 


These  equations  are  reduced  to  a  single  equation  in  p  and  the  resulting 
equation  is  solved  by  Newton’s  method  of  iteration. 


(4)  Numerical  Determination  at  lnt>  rior  Points 
Near  the  Boundary 

In  this  case  points  x,  and  x3  will  lie  between  Curve  I  and  the 
shock  w  ve,  but  point  x2  will  lie  outside  of  the  region  where  the  field 
quantities  are  to  be  determined. 

In  this  case  x'x  and  Xj,  together  with  the  field  quantities  at 
x\  and  Xj,  are  determined  by  the  same  formulas  as  for  the  general  interior 
point,  but  the  method  differs  for  determining  x'2  and  the  values  of  the 
field  quantities  at  x^. 


If  x I  as  computed  from 


[ ( v  ♦  c ) ( x 0 )  ♦  <  v  *  c) ( x  , )  ] 


A?  ♦  x 


is  greater  than  or  equal  to  x0,  the  whole  computat;  a  proceeds  a*  for 
the  general  interior  point,  with  the  exception  that 
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If  x'2  is  less  than  *0  ,  we  calculate  tf,  by  the  formula 


(  u  *  c)(x0)  ♦  ( t>  ♦  c)(xj) 

l  iTi 


t ( x)  -  v( xti) t(x0)  *  x0  -  X , 


(  v  *  c)(*0)  *  (  V  *  c)  (  X  J 


— .  (86) 


v{  XQ) 


In  this  case 


t\  -  tix0) 

At 


and  the  values  of  the  field  quantities  at  f2,  shown  in  Figure  15,  are 
found  by  using  Eq.  (60)  with  the  point  at  for  A,  and  *  for  B. 


FIG.  15  INTERIOR  POINT 

NEAR  THE  BOUNDARY 


The  f i ni te  di f ference  equations  for  determining  the  field 
quantities  at  x  are  then  given  by  Eqs.  (65)  through  (68)  together  with 
Eqs.  (69)  and  (70),  with  the  one  exception  being  that  A t  in  Eq.  (67) 
must  be  replaced  by  ( 1  -  \2)At. 

(5)  Numerical  I.etermi na t i on  at  an  Interior  Point 
Near  the  Shock  Wave 

This  is  the  final  type  of  point  which  must  be  considered  and 
is  shown  in  Figure  16. 
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FIG.  16  INTERIOR  POINT 

NEAR  THE  SHOCK  WAVE 


In  this  the  points  a'0  and  x  are  on  the  shock  Wjve,  but  the  points 
Xj  and  Xj  may  both  lie  to  the  right  of  the  shock  wave  or  only  x3  may  lie 
to  the  right.  In  any  event,  x2  will  lie  to  the  left. 

The  coordinate  of  x'2  is  found  by 


and 


if  x  Q  lies  to  the  ieft  of  x  or  x  j . 

If  xQ  lies  to  the  right  of  x0  or  x,  then  Xj  is  used  in  place  of 
x0  in  the  above  equat'on.  Again  the  field  quantities  at  x2  are  found  by 
Eq.  (60). 

If  x0  lies  to  the  right  ol  xJ(  then  x|  and  the  field  quantities 
at  xj  are  computed  as  for  the  general  interior  point.  If  x0  lies  to  the 
left  of  Xj,  then  x\  may  lie  either  to  the  left  of  or  right  of  x9. 

Thus  for  a'0  <  x^  xj  is  first  computed  from  the  formula 


( u(  x2  )  *  t'(  x0  )  ] 


57 


58 


The  values  of  the  field  quantities  at  tj  are  then  found  again 
from  Eq,  (60)  with  the  appropriate  but  obvious  substitutions. 

This  completei  the  method  of  computation  at  all  points  except 
the  initial  point. 

( 6 '  Initial  Values 

In  order  to  determine  the  initial  values  we  use  the  boundary 
conditions  and  the  jump  conditions  across  the  shock  wave,  p  x)  or  v(x) 
is  known  at  the  initial  point.  Then  Eqs.  (81),  (82),  and  (84)  can  be 
used  to  determine  the  remaining  three  of  the  four  quantities  U(x),  v(x), 
p(x),  and  p( x) .  Again  the  equations  are  reduced  to  a  single  equation  in 
p(x)  and  the  resulting  equation  is  solved  by  the  method  of  Newton.  This 
method  has  been  programmed  into  the  digital  computer  Burroughs  5500  and 
IBM  7090. 


Method  of  von  Neumann  and  Bichtmyer 


The  equations  of  motion  and  the  constitutive  relationships  can  be 
rewritten  as 

'<fcr  3 


P  V 


dx 


*x 


(-P  *  *.«> 


(87) 


P_ 

P 


dv 


(88) 


ftp) 


(89) 


2G 


D..  -  Qy 

2  GT 


3v 

Jx 


(90) 


These  equations  apply  only  when  lO^I  >  y.  In  these  equations,  the  dots 
ere  the  time  derivatives  along  the  particle  path.  A  simple  one- pa rameter 
relaxation  model  has  been  assumed,  for  purposes  of  discussion.  Other 
models  can  be  treated  in  tnis  way. 

One  of  the  methods  of  integration  of  the  partial  differential  equa¬ 
tions  of  the  type  expressed  in  Eqs.  (87)  to  (90)  is  the  method  of 
von  Neumann  and  Bichtmyer  \Befs.  14,15),  sometimes  called  the  “ Q"  method. 
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In  this  method,  finite  difference  forms  of  the  differential  equations  are 
used.  The  necessity  of  locating  the  shock  position  and  of  applying  the 
appropriate  jump  conditions  across  these  shock  waves  is  cl. initiated  by 
introducing  a  stress  term  defined  as  artificial  viscous  stress  Q.  This 
stress  is  proportional  to  the  gradient  of  the  particle  velocity  and  has 
the  effect  of  spreading  the  shock  discontinuities  over  a  few  cells  of  the 
finite  difference  network.  The  introduction  of  the  viscous  stress  intro¬ 
duces  an  error  in  addition  to  the  error  caused  by  the  finite  difference 
approximation  of  the  differential  equation.  However,  this  approximation, 
as  used  and  demonstrated,  is  quite  good  and  very  useful  in  solving  many 
practical  problems.  The  amount  of  error  can  be  controlled  by  controlling 
the  size  of  finite  difference  network  and  by  the  minimal  use  of  Q.  The 
use  of  Q  can  be  kept  to  a  minimum  by  introducing  the  artificial  viscosity 
stress  only  at  the  shock  fronts.  Even  though  the  shock  discontinuities 
are  spread  over  a  few  cells  of  a  finite  difference  network,  the  jump  con¬ 
ditions  across  the  shock  waves  are  still  satisfied  within  a  degree  of 
accuracy  that  can  be  controlled  by  the  size  of  the  finite  difference  net¬ 
work  and  the  magnitude  of  Q.  With  this  Q,  Eq.  (87)  changes  to 


pi 


3 

&  *  9 


(91) 


The  usual  expression  used  for  Q 


i  s 


1  dp 


where  A  is  a  constant.  Then 


(92) 


(93) 


It  is  the  stress  term  proporional  to  dv/dx  that  provides  dissipation 
and  spreading  of  the  shock'  wave. 


d.  Modified  Q  Method 

In  the  problem  of  wave  propagation  in  st ress- rel i*i ng  solids,  the 
dissipation  is  built  in  through  the  constitutive  relations.  This 


80 


dissipation,  which  comes  through  the  physical  mechanism  of  stress  relaxa 
tion,  can  take  the  place  of  the  artificially  introduced  Q  in  the  method 
of  von  Neumann.  This  can  be  seen  in  the  following  way.  Equation  (90^ 
can  be  simplified  to  yield 

D>%  -  2GT  ^  -  TDxt  *  By  .  (94) 

This  equation  and  Eq.  (87)  yield 

d  /  '  dv  \ 

pv  =  -  -p  *  By  -  TD  *  2GT  —  .  (95) 

0*  \  Ox j 

Equation  (95)  contains  a  term  proporti oral  to  3t’/3x,  just  as  in  Eq.  (93) 
This  can  provide  the  spreading  of  the  plastic  shock  wave.  Thus  we  have 
a  dissipation  term  from  the  physical  mechanism  rather  than  an  artificial 
one.  However,  in  starting  the  process  of  integration  a.i  artificial  Q  is 
needed  until  the  calculated  stress  deiiator  terms  are  sufficient  to  take 
care  of  stabilising  the  procedure  of  numei i, integration. 

The  stability  and  convergence  of  such  methods  are  discussed  in 
Refs.  15  and  16. 

Finite  difference  equitions  and  the  coding  used  in  numerical  inte¬ 
gration  are  similar  to  those  discussed  by  Wilkins  (Ref.  17). 


11 


4,  Cecay  of  the  Elastic  Precursor 


The  decay  of  pressure  at  the  elastic  shock  front  can  be  calculated 
in  a  stress- relaxing  material  under  restrictive  assumptions.  The  model 
for  this  discussion  will  be  the  one-dimensional  Zener  t vpe  of  macroscopic 
relaxing  material.  As  discussed  in  Section  III-2-bf 11  ,  the  constitutive 
relation  for  this  material  will  be 


P,  "  r-2P 


F'lp,ax) 


A  restrictive  assumption  which  makes  some  analytical  calculation 
possible  is  that  the  px  -  p  relationship  in  the  unrelaxed  (t.e.,  elastic) 
material  is  given  by 


P, 


*  A 


(96) 


This  relationship  is  similar  to  the  hydrodynamic  type  of  relation¬ 
ship.  It  will  introduce  a  very  serioi-s  restriction  on  the  application  of 
the  results;  however,  a  trend  might  be  indicated. 

Along  the  leading  characteristic  a, 


—x 

dt 


♦  pc 


dv 

Tt 


-F' 


(97) 


and  the  path  of  the  characteristic  OL  is  defined  by 


dx 

dt 


v 


♦  c 


(98) 


The  constant  A  can  be  written  as 

A  -  p0t/02  (99) 

where  (/,  is  a  quantity  with  dimensions  of  velocity. 

At  the  front- running  shock  wave  in  a  relaxing  material  the  jump 
conditions  are  elastic.  Then  the  front  .unning  shock  velocity  V  can  be 
calculated  from  the' jump  condi tio  .s. 
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r2 


P . 


■(>  '  7) 


(100) 


From  Eqs.  (98)  and  (99) 


"02 


Hence 


P, 


4-  7) 


(101) 


u  =  un 


(102) 


Then  the  jump  condition  expressing  the  conservation  of  mass  yields  the 
Darticle  velocity 


4  -  7) 


(103) 


dP,  (  Pouo 


dp  \  p 


I  •  ■  • 


.  (104) 


The  path  of  the  leading  characteristic  ot  can  be  calculated  from  Eq.  (97) 


'dx 

LGt 


V  *  c 


Un 


(105) 


Thus  the  leading  characteristic  ®  and  the  front- running  (or  elastic)- 
shock  wave  coincide.  Along  the  elastic  shock  wave  the  compatibility 
equation  [Eq.  (97)]  corresponding  to  the  at  characteristic  is  valid.  By 
usi'if  the  expression  for  c  from  Eq.  (104),  Eq.  (97)  becomes 


dv 

77  +  p  ^  Tt 


(106) 


Along  the  front- running  shock  wave  the  jump  conditions  yield 


F,  B  P%Vv  =  psU0v 


(107) 
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Then  t'.q.  (106)  becomes 


dP, 

dt 

It 

( 1  OR ) 

This  equation 

gives  the 

rate  of  decay  at  the  wav'-  front. 

If  F'  is  given 

F‘ 

P,  ~  P. 

T 

(  109) 

'Then 

dP , 

dt 

P,  *  P> 

2  T 

(110) 

dP, 

-  -  1* 

2  T 

(111) 

Px  -  P, 

and 

P> 

-  Pi  ♦  Be~’/2T 

(112) 

where  B  is  an  arbitrary  constant.  The  constants  B  and  T  can  be  deter¬ 
mined  from  experimental  results.  If  F'  is  given  by 


and 


then 


and 


F'<p,at) 


P,  ~ 

’  Ti 


(113) 


dpx 

dt 


Pi 


p  *  [p  t  ( 0 )  -  p]e 


- «/  r. 


(114) 

(115) 


dP, 

~dt 


_  -  t  /  T  _ 

[p^O)  -  p)  e  2  -  p)  .  (116) 
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I  .  €  .  , 


111 

d  t 


(117) 


1  P,  [p ,  f  0)  -  pi 


-t/T 


2  r, 


2  T, 


2  , 

2  r, 


The  solution  to  this  differential  equation  is 


Pm 


If  the  initial  amplitude  of  the  precursor  is  different  from  the  initial 
yield  value  p^( 0),  then  A (  is  not  zero  and  the  relaxation  constant  T j 
tekes  care  of  the  time  required  to  relax  toward  the  yield  value  and  the 
formation  of  a  two-wave  structure.  This  two-wave  structure  is  continually 
changing  with  decreasing  amplitude  of  the  precursor  toward  the  st  tic 
value.  This  is  taken  care  of  by  the  relaxation  constant  Tj. 

These  calculations  are  similar  to  the  calculation  by  Taylor,  and  any 
evaluation  of  relaxation  phenomena  in  this  way  is  baoed  on  some  very  re¬ 
strictive  assumptions,  as  discussed  earlier. 

-J*- 

Any  accurate  evaluation  of  the  relaxation  phenomena  must  be  based  on 
realistic  p x  -  p  relationships  and  the  models  based  on  continuum  mechanics. 
At  present  this  can  be  done  only  by  numerical  techniques.  Furthermore, 
these  methods  provide  some  freedom  in  the  choice  of  parameters  in  the  re¬ 
laxation  model,  and  thus  the  experimental  results  including  the  flow  behind 
the  wave  f^ont  can  be  interpreted  more  accurately. 


5.  Discussion  of  Some  Numerical  Results 


a.  Wave  Calculations 

Experimental  observations  in  hard  and  soft  aluminum  indicate  a  two- 
wave  structure  at  very  early  times;  the  amplitude  of  the  elastic  wave 
decreases  with  time.  The  two-parameter  macroscopic  relLxation  model  was 
used  in  a  modified  Q  code  to  correlate  the  experimental  results.  The 
properties  of  aluminum  used  in  the  calculation  are  as  follows; 

P0  *  2.785  gm/cc 

p  =  a(p/pc  -  1)  ♦  Mp/p0  "  l)2  ♦  c(p/pn  -  l)3 
a  =  0.755  Mbar 
b  =  1.29  Mbar 
c  =  1.197  Mbar 
7  =  0.0025  Mbar 
G  =  0.287  Mbar 

The  relaxation  constants  7\  and  T2  were  introduced  earlier  in 
Section  III-2-b(l ) .  T}  controls  the  rate  of  formation  of  e  two-wave 
structure;  expeiimental  observation  indicates  that  it  should  be  very 
small.  T2 . governs  the  rate  of  decay  of  the  precursor  amplitude. 

Th *  values  of  T j  and  T 2  used  in  the  calculation  are  0.005  and  3.72 
psec,  respectively.  These  values  were  efiosen  by  trial  and  error  to 
match  approximately  the  observed  precursor  decay — they  do  not  represent 
a  best  fit.  Figure  17  shows  the  decay  of  the  precursor  amplitude  with 
the  penetration  depth.  A  comparison  of  observed  and  computed  precursor 
decay  is  shown  in  Fig.  18. 

The  curve  in  Fig.  19  gives  the  variation  of  crg  with  distance  for 
selected  times  when  the  boundary  is  subjected  to  a  fixed  velocity  rather 
than  a  fixed  stress.  This  curve  was  computed  by  the  characteristic  method 
for  a  one-parameter  model;  hence  the  profiles  do  not  appear  as  two  separate 
waves  as  observed  experimentally.  Note,  however,  the  pronounced  stress 
relaxation  at  the  boundary.  Computations  like  this  are  needed  {or  the  two- 
parameter  model,  so  that  a  more  direct  comparison  with  experiment  can  be 
made. 
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FIRST  WAVE  AMPLITUDE 
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IG.  19  STRESS  PROFILES  FOR  CONSTANT  VELOCITY  INPUT  AT  THE  INTERIOR  BOUNDARY 

b.  Computation  from  the  Characteristic  Code 

In  the  problem  of  stress  relaxation  we  have  a  dissipation  built  in 
through  the  constitutive  relations.  In  the  modified  0  method  described 
in  Section  III-3-d,  the  artificial  dissipation  has  been  kept  to  a  minimum; 
it  is  confined  to  the  starting  of  the  process  of  computation.  However, 
one  weakness  of  the  method  is  that  the  elastic  shock  front  is  not  clearly 
defined  and  the  elastic  wave  amplitude  cannot  be  predicted  exactly.  To 
overcome  these  deficiencies,  the  characteristic  code  is  the  most  accurate 
method  available.  The  error  in  this  method  is  due  only  to  the  approximation 
of  the  differential  equation  by  a  finite  difference  scheme.  This  can  be 
controlled  by  suitably  selecting  the  finite  difference  mesh  work,  depend¬ 
ing  on  the  desired  degree  of  accuracy. 


Early  numerical  calculations  were  done  for  quartz  for  purposes  of 
illustration  of  the  stress- relaxation  properties.  (At  the  time  these 
calculations  were  made,  experimental  results  on  aluminum  were  not  avail¬ 
able.)  A  one-parameter  macroscopic  relaxation  model  was  used. 

The  p  ~  p  relationship  for  quartz  was  assumed  to  b£ : 

p  =  0.3919  (p/pQ  -  1)  Mbar 

P0  ~  2.5  gm/cc 


Other  parameters  used  were: 

y  =  0.002  Mbar 

G  =  0.37  Mbar 

T  -  0.  03  naec. . 

Figure  20  is  an  x  -  t  diagram  showing  the  positions  of  the  left 
boundary  and  the  position  of  elastic  shock  front  with  time.  Figure  21 
illustrates  the  distribution  of  pz  with  the  Eulerian  coordinate  x  for 
different  values  of  time  t.  In  the  computation  the  stress  on  the  left 
boundary  was  kept  constant  at  0.1  Mbar.  The  curves  indicate  the  results 
obtained  by  the  method  of  characteristics;  the  points  illustrate  the 
values  obtained  by  the  modified  Q  method. 

Figure  22  illustrates  the  decay  of  the  elastic  shock.  Figures  23 
and  24  illustrate  the  distributions  of  the  particle  velocity  and  the 
specific  volume  with  the  Eulerian  coordinate  x  for  different  values  of 
time  t. 
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FIG.  20 
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FIG.  22  DECAY  OF  PRECURSOR  STRESS  WITH  PROPAGATION  DISTANCE 
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FIG.  23  PARTICLE  VELOCITY  PROFILES  AT  DIFFERENT  TIMES 


FIG.  24  SPECIFIC  VOLUME  PR0FILE5  AT  DIFFERENT  TIMES 
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Dislocation  Model 


Figures  25-29  illustrate  the  results  obtained  from  numerical  calcula 
tions  (von  Neumann  and  Richtmyer  method)  when  the  dislocation  model  for 
stress  rc'axation  was  used  [see  Section  1 1 1 -2-b- (2 ) ] •  Again,  the  cal¬ 
culations  were  performed  only  for  purposes  of  illustration.  Iron,  which 
exhibits  pronounced  stress  relaxation,  was  used  for  calculations.  The 
following  values  were  used  for  the  constants  in  the  constitutive  relation. 
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FIG.  25  SHOCK  PROFILES  SHOWING  THE  EFFECT  OF  DISLOCATION  MULTIPLICATION 
Profiles  in  iron  at  0.25  and  0.45  pitc. 
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FIG.  26  SHOCK  PROFILES  SHOWING  THE  EFFECT  OF  INCREASING  g. 
Profiles  in  iron  at  0.25  and  0.45  /xsec. 
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FIG.  27  SHOCK  PROFILES  SHOWING  THE  EFFECT  OF  INCREASING  THE  INITIAL 
DENSITY  OF  PINNED  DISLOCATIONS.  Profile*  in  iron  ot  0.25  ond  0.45 usee. 
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FIG.  2B  PRECURSOR  DECAY  FOR  VARIOUS  VALUES  OF  B, 
THE  CONSTANT  IN  THE  GiLMAN-JOHNSTON  MODEL 


FIG.  '  PRECURSOR  DECAY  FOR  VARIOUS  VALUES  OF  THE  INITIAL 
\BILE  DISLOCATION  DENSITY 
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Figure  25  shows  the  d  i  st.  rihut  ion  of  with  the  Kulerian  coordinate  x  for 

two  specific  times:  t  r  0.25  /i sec  and  t  -  0.45  psec.  This  figure  illus- 
0 

trates  the  effect  of  dislocation  multiplication;  the  open  circles  are  for 
no  dislocation  multiplication,  and  the  dots  are  for  g  =  6000  (Mbar-ptsec)  1 
and  bN t 0  =  7.5  cm*1.  The  effects  of  dislocation  multiplication  are  a  more 
pronounced  elastic  precursor  and  a  steeper  plastic  shock  wave. 

Figure  26  illustrates  the  effect  of  changing  g.  Two  values  of  g 
were  used:  g  -  60CO  (Mbar-ptsec)  1  and  g  =  20,000  (Mbar-pisec)  *.  A  value 
of  25  cm-1  was  used  for  6A't0  in  both  cases.  The  higher  value  of  g  reduced 
the  precursor  amplitude  significantly. 

Figure  27  illustrates  the  variation  of  bN  t  Q  for  a  fixed  value  of  g 
of  6000  (Mbar-pisec)  '.  Two  values  of  bN (  0  , 

bNt0  =  7.5  cm"1 

and 

bNt0  -  25.0  cm" 1  are  compa  red . 

The  larger  value  of  has  the  lower  amplitude  of  elastic  precursor. 

Figures  28  and  29  illustrate  the  decay  of  the  amplitude  of  the 
elastic  shock  wave  with  time  for  different  values  of  B  and  N  n. 

a  0 

These  parameter  studies  illustrate  the  effect  of  various  material 
constants  on  shock  propagation. 

d .  Attenuation  of  Peak  Pressure 

One  of  the  important  problems  of  practical  interest  is  th.  attenua¬ 
tion  of  peek  pressure  of  a  pulse  of  finite  duration.  A  study  of  this 
phenomenon  requires  a  knowledge  of  the  unloading  behavior  and  the 
constitutive  relations  during  unloading.  The  constitutive  relations  for 
'oaiitng  are  now  understood  in  some  cases  where  stress  relaxation  is  pre¬ 
sent,  however,  there  are  no  experimental  data  available  on  unloading 
under  conditions  such  that  stress  relaxation  is  important. 


To  obtain  some  information  regarding  the  trend  that  the  decay  of  peak 
pressure  might  follow  when  stress  relaxation  is  present,  the  following 
calculation  has  been  performed.  It  has  been  assumed  that  the  relaxing 
constitutive  relation  is  the  same  in  loading  as  well  as  in  unloading. 

The  material  and  the  constants  used  are  the  same  as  in  the  case  of  the 
computation  with  the  characteristic  code.  A  pressure  of  0.1  Mbar  is 
applied  at  time  t  -  0;  the  subsequent  loading  history  is  shown  in  Fig.  30. 
This  figure  illustrates  the  decay  of  peak  pressure  with  time.  The  dashed 
line  co.responds  to  e  las t i c - p last ic  theory  without  stress  relaxation. 


FIG.  30  CALCULATED  DECAY  OF  PEAK  PRESSURE  IN  ALUMINUM 
WITH  AND  WITHOUT  STRESS  RELAXATION 


6.  Summary 


In  summary,  the  major  results  of  the  theoretical  part  of  this  con¬ 
tract  are  the  devel npmcn t  of  a  mathematical  model  of  stress  relaxation, 
together  with  numerical  results  on  the  rate  of  decay  of  the  elastic  pre¬ 
cursor  wave,  and  some  computed  wave  profiles.  The  numerical  results  have 
been  correlated  with  experimental  results  and  it  is  shown  that  the  two- 
parameter  model  of  stress  relaxation  predicts  experimental  results  more 
accurately  then  the  one  parameter  model.  The  computations  have  been  per¬ 
formed  using  the  method  of  characteristics  and  the  so-called  “Q"  method 
of  von  Neumann  and  Richtmyer.  The  reasons  for  the  expenditure  of  effort 
on  the  method  of  characteristics  are: 

1.  The  method  is  r.ora  accurate  thin  the  “Q”  method.  It  was  there¬ 
fore  possible  to  study  the  effect  of  stress  relaxation  in  the  absence  of 
artificial  viscosity,  then  to  use  the  “Q”  method  to  study  stress  relaxation, 
and  to  compare  the  difference  in  the  results  obtained  by  the  two  methods. 

2.  The  “Q"  method  is  highly  inaccurate  for  situations  involving  re¬ 
flected  waves;  however,  results  obtained  by  the  method  of  characteristics 
can  be  used  to  indicate  correction  terms  which  could  be  added  to  “Q” 
method  codes. 

Another  result  of  the  present  work  is  that  the  term  involving  stress 
relaxation  provides  a  real  viscosity  and  hence  the  artificial  viscosity 
term  may  be  omitted  from  “Q"  codes  when  stress  relaxation  is  present.  It 
should  bo  pointed  out  that  the  stress  relaxation  models  employed  here  are 
most  applicable  to  metals  and  it  is  quite  possible  that  thermodynamic 
effects  will  be  of  greater  importance  in  plastics  than  stress  relaxation; 
this  point  requires  further  study. 

It  is  to  be  noted  that  the  basic  differential  equations  used  in  the 
present  work  difer  from  those  used  in  the  PUFF  codes  only  in  the  con¬ 
stitutive  relations  which  are  used.  It  is  not  a  difficult  matter  to  in¬ 
corporate  the  stress  relaxation  model  into  a  PUFF  code.  Some  work  on 
expansion  of  the  models  and  more  correlation  with  experimental  results 
should  be  accomplished  before  these  models  are  used. 
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APPENDIX  I 

SHOTS  FIRED  DURING  THIS  PROGRAM 


All  the  shots  fired  under  the  present  contract  are  listed 
with  an  indication  of  the  puipose  for  firing  each  shot .  Some 
were  of  a  development  nature  and  did  not  contribute  significan 
points. 


below, 
shot  s 
t  data 


SHOT  NO. 


TYPE  Of 
GAGE 


APPENDIX  I 


PURPOSE 


12,810 

Quart i 

Calibration  of  quart*  gages 

12,812 

Quart* 

Calibration  of  quart*  gage* 

1 2 , 8 1  3 

Quart* 

Calibration  of  quart*  gages 

12,811 

2  Quartz 

Development  of  technique  to  mount  ouarti  in  projectile 

12,053 

Quartz 

Look  for  elaatic  rave  in  Plexiglas 

12,052 

Quart* 

Look  at  elastic  wave  in  stock  202  1  aluminum 

12,050 

Mangsnin  ^  *  * 

Look  at  elaatic  Bare  in  stock  2024  aluminum 

12,051 

u  -(b) 

Mangani n 

Look  at  elastic  wave  in  stock  2024  aluminum 

12,172 

Quart* 

Study  elaatic  wave  in  annealed  2024  aluminum 

12,168 

Quirt  * 

Study  elaatic  wave  in  annealed  2024  aluminum 

12,170 

Quartz 

Study  elaatic  wave  in  a,  sealed  2024  aluminum 

12,175 

Quart* 

Study  elaatic  wave  in  annealed  2024  tlummum 

12. 160 

Quart* 

Study  elastic  wave  in  annealed  2024  aluminum 

12.807 

Quart* 

Study  elaatic  wave  in  hardened  2024  aluminum  and  look  for  plastic 

wave  decay 

12,808 

Quartz 

Studv  elaatic  wave  in  hardened  2024  aluminum  and  look  for  plastic 

wave  decay 

12,173 

Quart* 

Study  elaatic  wave  in  hardened  2024  aluminum 

12,171 

Quart* 

Study  elaatic  wave  in  hardened  2024  aluminum 

12,17* 

Quirt* 

Study  elastic  wave  in  hardened  2024  aluminum 

12.311 

Manganin 

Study  wave  profile  in  polyethylene 

12,512 

2  Maogani 

Study  save  profile  in  polyethylene 

12,513 

2  Manganin 

Study  wave  profile  in  polyethylene 

12,51* 

2  Manganin 

Study  weve  profile  m  polyethylene 

12.515 

2  Manganin 

Study  wave  profile  in  polyethylene 

12.516 

2  Manganin 

Study  wave  profile  in  polyethylene 

12, 517 

1  Manganin 

Study  wave  profile  in  polyethvlene 

12.80* 

2  Manganin 

Study  wa«e  profile  in  polyethylene 

12.805 

2  Manganin 

Study  wave  profile  in  polyethylene 

12.806 

2  Manganin 

Study  save  profile  in  polyethylene 

Hh|mi a  tirt  fr**  AI  plai*  ky  OJ(01-i»ck  Mylar  foil  ^ottoA  »•  C-?  apoay. 

lu|Hii  oi  ro  silk  0. 061  -  tor  K  Mylar  »4  Mrrourfif  bv  atuuiaua  block*. 
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APPENDIX  II 


CHARACTERISTIC  EQUATION  BASED  ON  THE  EXISTENCE 
THEOREM  OF  CAUCHY -KOMAlEfSKI 


APPENDIX  II 


CHARACTERISTIC  EQUATION  BASED  ON  THE  EXISTENCE 
THEOREM  OF  CAUCHY -KOIALEISK I  (Ref.  18) 


We  begin  with  the  equations  expressing  conservation  of  mass  and 
momentum 


dcr 
B x 


a) 
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(2) 


We  write  the  constitutive  relations  in  the  form 
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These  equations  may  be  conveniently  rewritten  in  the  form 
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2  dv 
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|_2G 


3c  1  SSj 

3  9*  1 2 G  3f  3* 


=  g(5,  -  0Y) 


(8) 


Altiiough  many  equivalent  definitions  of  characteristics  are  possible, 
we  adapt  here  the  one  based  on  the  existence  proof  of  Cauchy-Kowa lewski . 
The  reason  for  this  choice  is  that  it  gets  us  to  the  desired  goal  along 
the  most  direct  path. 

We  omit  obvious  details  and  proceed  to  look  for  curves  f(x,t)  -  con¬ 
stant,  along  which  the  set  of  Kqs.  (5)  -  (8)  cannot  be  solved  as  an 
algebraic  system  for  the  partial  derivatives  3Sj/3/,  3 p/of,  9v/3 /  and 

Si  nc  e 

*1  'Ifi 

9x  3/  9x  3 f  9 1 

etc.,  we  obtain  the  following  set  of  equations: 

3S, 


3/  3 p  9/  3v  /3 f  *  3 /  \ 

3x  3 f  3x  9/  \3f  3 z  J 

3A  /3/  3/  \  9i>  3/ 

9 f  \3 1  3 xj  9 f  3* 

3p  /9/  3/  \  dj  ( 3/ 

9/  \  3(  3.t  y 


2  H  3/ 

3  '7  3, 


1  pi 

2c;  7  \3f 


9/  /  3/  ^  3/  \ 

dp  3 f  \<lt  9*  / 


3f  9f\ 

♦  e  —  -  guS.  -  GK)  =  0 

9x  )  1 


The  matrix  of  the  coefficients  of 


AS 


3y 

,  *“  ,  “ *•  ,  and 

3 1  9/  9/  9/ 


80 


.tte&as  A  ..  .  -  it  ■ 


and  the  determinant  of  this  matrix  is 


3/Y 

-  I 

f  3J. 

,2fY  -*( 

pA2  „  £ 

df_ 

(IlY 

\3t 

ox) 

2  G 

[3, 

dx)  3  1 

\^x)  2  G 

dp 

\^x) 

There  will  be  no  solution  lor  etc.  if  this  determinant  vanishes. 

On  /  =  constant  d//3f  =  - (<3//3x ) (dx/dt  )  so  that  the  determinant  vanishes 
when 


dx 

7t 


(9) 


The  four  roots  given  by  l£q .  (9)  are  the  characteristics  of  the  govern¬ 
ing  differential  equations. 
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APPENDIX  111 


COMPATIBILITY  C8MITI0KS 


The  governing  differential  equations  are: 


(1) 


(2) 


3s,  35 


l  AG  3t> 


3t  3*  3 


—  -  2Cg(5j  -  6Y) 


3  p  3  p 

17*  <'* 


\  3 1  dx) 


(3) 


(4) 


The  four  characteristics  are: 


dx 


where 


77  " 


df  AG 

dp  V 


If  we  multiply  Eq.  (I)  by  c,  add  Er(.  (3)  and  subtract  Eq.  (4),  we 
obt  a i n 


3  3 

—  (5,  -p)  ♦  <i  ♦  c)  — 
3|  • 


/3i  3u\  AC  3v  3v 

(S,  •  .  -j-  T  -  -Pi  </•>  jj  ■  2T«(S,  -  H»-> 


K 


Then  using  Eq.  (2)  we  get 


9  3 

~  (S.-p)  +  (v+c)  —  (S  -  p ) 
Of  ox  1 


-  pc  I  ~  4  (  V  +  c  ) 


2Gg(5,  -  6 >Y) 


Thus  we  have 


a  long 


a  1  ong 


d(S!  -  p)  -  pcdv  -  2Gg (S j  -  6Y)dt  =  0 


d{S  j  -  p)  >  prdv  -  2Gg(Sj  -  ct')dt  =  0 


By  using  Eq .  (2)  in  Eq.  (3)  we  get 


T'l 

dS  t  ♦  —  dp  -  2Cg(S,  -  r))dt 


along 


and  Eq.  (i)  can  be  wvitten  directly 


1 1  v  as 


dp  -  c ^dp  ~  0 


•  1  ong 


Equations  (5),  <  b ) .  (7).  and  (8)  a t  the  compatibility  equations 
lor  a  stress-relaxing  medium. 
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